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ABSTRACT 

X 

A pseudo-Nambu-Goldstone boson, with a potential of the form V(4>) = A 4 [l ± cos((j)/ /)], 
can naturally give rise to an epoch of inflation in the early universe, if / ~ Mpi and A ~ Mqut- 
Such mass scales arise in particle physics models with a gauge group that becomes strongly in- 
teracting at the GUT scale. We explore the particle physics basis for these models, focusing 
on technicolor and superstring theories, and work out a specific example based on the multiple 
gaugino condensation scenario in string/supergravity theory. We study the cosmological evolution 
of and constraints upon these models numerically and analytically. To obtain a sufficiently high 
post-inflation reheat temperature for baryosynthesis to occur we require / > 0.3M p /. The primor- 
dial density fluctuation spectrum generated by quantum fluctuations in <ft is a non-scale-invariant 
power law, P(k) tx k Us , with n s ~ 1 — (M|,^/87r/ 2 ), leading to more power on large length scales 
than the n s = 1 Harrison-Zeldovich spectrum. We pay special attention to the prospects of us- 
ing the enhanced power to explain the otherwise puzzling large-scale clustering of galaxies and 
clusters and their flows. We find that the standard cold dark matter model with < n s < 0.6 
could in principle explain this data. However, the microwave background anisotropies recently 
detected by COBE imply such low primordial amplitudes (that is, bias factors b$ > 2) for these 
CDM models that galaxy formation would occur too late to be viable and the large-scale galaxy 
flows would be too small; when combined with COBE, these each lead to the constraint n s > 0.6, 
hence / > 0.3Mpi, comparable to the bound from baryogenesis. For other inflation models which 
give rise to initial fluctuation spectra that are power laws through the 3 decades in wavelength 
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probed by large scale observations, such as extended inflation and inflation with exponential po- 
tentials, our constraint on n s is tighter, n s > 0.7. Combined with other constraints (which imply 
n s < 0.77 — 0.84), this leaves little room for most extended inflation models. Chaotic inflation 
models with power law potentials have n s > 0.95 through this band and so are not affected. 
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I. INTRODUCTION 



In recent years, the inflationary universe has been in a state of theoretical limbo: it is a 
beautiful idea in search of a compelling model. The idea is remarkably elegant [1]: if the early 
universe undergoes an epoch of quasi-exponential expansion during which the Robertson- Walker 
scale factor a(t) increases by a factor of at least e 60 , then a small causally connected region grows 
to a sufficiently large size to explain the observed homogeneity and isotropy of the universe, to 
dilute any overdensity of magnetic monopoles or other unwanted relics, and to flatten the spatial 
hypersurfaces, Q = 87rGp/3H 2 — > 1. As a bonus, quantum fluctuations during inflation can 
causally generate large-scale density fluctuations, which are required for galaxy formation[2]. 

During the inflationary epoch, the energy density of the universe is dominated by the (nearly 
constant) potential energy density V{4>) associated with a slowly rolling scalar field (f>, the in- 
flaton[3]. To satisfy cosmic microwave background radiation (CMBR) anisotropy limits on the 
generation of density fluctuations, the potential of the inflaton must be very flat. Consequently, 
the field (f> must be extremely weakly self-coupled, with effective quartic self-coupling constant 
satisfying A^ < 10~ 12 — 10~ 14 in most models [4]. 

Density fluctuations in inflation are thus a blessing for astronomers but a curse for particle 
physicists, because the theory must contain a very small dimensionless number. Attitudes con- 
cerning this problem vary widely among inflation theorists: to some this represents unacceptable 
'fine tuning'; to others, it is not an issue of great concern because we know there exist other 
small numbers in physics, like lepton and quark Yukawa couplings gy ~ 10 -5 and the ratio 
M wea k/Mpi ~ 10~ 17 . Partly as a consequence of the latter view, in recent years, it has become 
customary to decouple the inflaton completely from particle physics models, to specify an 'inflaton 
sector' with the requisite properties, with little or no regard for its physical origin. 

Nevertheless, it is meaningful and important to ask whether such a small value for is 
in principle unnatural. Clearly, the answer depends on the particle physics model within which 
4> is embedded and on one's interpretation of naturalness. A small parameter A is said to be 
"technically natural" if it is protected against large radiative corrections by a symmetry, i.e., if 
setting A — > increases the symmetry of the system [5]. For example, in this way, low energy su- 
persymmetry might protect the small ratio M wea k/Mpi. However, in technically natural inflation 
models, the small coupling A^, while stable against radiative corrections, is itself unexplained, 
and is generally postulated (i.e., put in by hand) solely in order to generate successful inflation. 
Technical naturalness is a useful concept for low energy effective Lagrangians, like the electroweak 
theory and its supersymmetric extensions, but it points to a more fundamental level of theory for 
its origin. Since inflation takes place relatively close to the Planck scale, it would be preferable to 
find the inflaton in particle physics models which are "strongly natural", that is, which have no 
small numbers in the fundamental Lagrangian. 

In a strongly natural gauge theory, all small dimensionless parameters ultimately arise dy- 
namically, e.g., from renormalization group (or instanton) factors like exp( — 1/a), where a is a 
gauge coupling. In particular, in an asymptotically free theory, the scale Mi, at which a loga- 
rithmically running coupling constant becomes unity, is small, Mi ~ M2e -1 / a , where M2 is the 
fundamental mass scale in the theory. In some models, the inflaton coupling A^ arises from a 
ratio of mass scales, A^ ~ (Mi/M 2 ) ra ; for example, in the models to be discussed below, n = 4. 
As a result, in such models, A^ is naturally exponentially suppressed, A^ ~ e~ n l a . 

An example of this kind, namely, a scalar field with naturally small self-coupling, is provided 
by the axion [6], a light pseudoscalar which arises in models introduced to solve the strong CP 
problem. In axion models, a global U(l) symmetry is spontaneously broken at some large mass 
scale /, through the vacuum expectation value of a complex scalar field, (<£) = / exp(ia/f). 
(In this case, $ has the familiar Mexican-hat potential, and the vacuum is a circle of radius /.) 
At energies below the scale /, the only relevant degree of freedom is the massless axion field 
a, the angular Nambu-Goldstone mode around the bottom of the $ potential. However, at a 
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much lower energy scale, the symmetry is explicitly broken by loop corrections. For example, 
the QCD axion obtains a mass from non-perturbative gluon configurations (instantons) through 
the chiral anomaly. When QCD becomes strong at the scale Aqcd ~ 100 MeV, instanton effects 
give rise to a periodic potential of height ~ Aq CD f° r the axion. In 'invisible' axion models [7] 
with canonical Peccei-Quinn scale fpQ ~ 10 12 GeV, the resulting axion self-coupling is extremely 
small: A a ~ (Aqcd/ /pq) 4 ~ 10~ 52 . This small number simply reflects the hierarchy between the 
QCD and Peccei-Quinn scales, which arises from the slow logarithmic running of oqcd- 

Pseudo-Nambu-Goldstone bosons (PNGBs) like the axion are ubiquitous in particle physics 
models: they arise whenever an approximate global symmetry is spontaneously broken. We there- 
fore choose them as our candidate for the inflaton: we assume a global symmetry is spontaneously 
broken at a scale /, with soft explicit symmetry breaking at a lower scale A; these two scales com- 
pletely characterize the model and will be specified by the requirements of successful inflation, 
namely, a sufficent number of e-folds of inflation, sufficient reheating, and an acceptable amplitude 
and spectrum of density fluctuations. The resulting PNGB potential is generally of the form 

y(0) = A 4 [l±cos(iV0//)]. (1-1) 

We will take the positive sign in Eq.(l.l) (this choice has no effect on our results) and, unless 
otherwise noted, assume N = 1, so the potential, of height 2A 4 , has a unique minimum at 4> = irf 
(we assume the periodicity of 4> is 27r/). In a previous paper [8] (hereafter Paper I), three of us 
showed that, for / ~ M Pl ~ 10 19 GeV and A ~ M GUT ~ 10 15 GeV, the PNGB field </> can drive 
inflation; in this case, the effective quartic coupling is ~ (A//) 4 ~ 10~ 13 , as required. In this 
paper, we study this class of models and their implications in greater depth. 

We note that, in some cases, the potential of Eq.(l.l) is the lowest order approximation to 
a more complicated expression. For inflation, the important ingredients are the height (~ A 4 ) 
and width (~ /) of the potential, and the curvature in the vicinity of its extrema, which is 
determined by m 2 = A 2 //. Thus, while our treatment will focus on the specific form (1.1), our 
conclusions hold for more general forms of the PNGB potential which have the same overall shape 
(that is, same height, width, and curvature at the extrema; in addition, we assume V((p) varies 
monotonically between (f> = and irf, that is, we ignore higher order ripples, which might affect 
the perturbation spectrum over a small range of wavelengths). 

In section II, we discuss the PNGB inflation scenario in the context of particle physics models. 
As noted above, a successful inflation scenario does not consist simply of a scalar field potential 
that does the trick; in addition, the parameters of the potential, in this case the requisite mass 
scales / and A, must have a natural origin in plausible particle physics models. PNGB potentials 
with these mass scales do arise naturally in particle physics models. For example, in the hidden 
sector of superstring (supergravity) theories, if a non-Abelian subgroup(s) remains unbroken, the 
running gauge coupling can become strong at the scale ~ 10 14 — 10 15 GeV; indeed, it is hoped that 
the resulting gaugino condensation may play a role in determining the string coupling constant 
and possibly in breaking supersymmetry [9]. (We note that, in such models, the only fundamental 
scale is the Planck scale, / ~ Mpi, and the lower scale A is generated dynamically.) In this case, as 
discussed in Section II, the role of the PNGB inflaton could be played by the "model-independent 
axion" (the imaginary part of the dilaton) [10]. 

In Secton III, we provide a detailed analysis of the cosmological evolution of the PNGB 
inflaton field. By and large, the numerical results therein confirm the analytic treatment of paper 
I. In addition, we also discuss in detail constraints on the mass scales arising from the requirement 
of sufficient reheating, the density fluctuation amplitude, and the requirement that inflation be 
probable in the sense of initial (and final) conditions. We also discuss the issue of initial spatial 
gradients in the inflaton field and how they may be damped out prior to inflation. 

In the standard lore of inflation, the adiabatic density fluctuations generated have a nearly 
scale-invariant Harrison-Zeldovich spectrum. This general statement can be violated, and an ar- 
bitrary perturbation spectrum 'designed', but at the cost of fine-tuning parameters of the inflaton 
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potential (or adjusting coupling constants in models with multiple scalar fields) [11]. One can 
imagine corrections to the PNGB potential which would allow this behavior, but in this paper we 
consider only the simplest model given by Eq.(l.l); in this case, we have no freedom to introduce 
features into the perturbation spectrum. Nevertheless, as discussed in section IV, in this model 
the fluctuations can deviate significantly from a scale- invariant spectrum: for / < 3M p j/4, the 

perturbation amplitude at horizon-crossing grows with mass scale M as (5p/p)hor ~ M m p l ^ 4Sn ^ . 
Thus, the primordial power spectrum for density fluctuations (at fixed time) is a power law, 
(\5p(k)/p\ 2 } ~ k ns , with spectral index n s ~ 1 — {M^JSnf 2 ). The extra power on large scales 
(compared to the scale-invariant n s = 1 spectrum) can have important implications for large-scale 
structure, of particular interest since the scale-invariant spectrum with cold dark matter (CDM) 
appears to have 8too little power on large scales. Other inflation models can also give rise to non- 
scale-invariant power law spectra. Therefore, in section IV, we discuss tests of non-scale-invariant 
power law initial spectra with adiabatic perturbations and CDM, including the galaxy angular 
correlation function inferred from deep photometric surveys, the CMBR anisotropy detected by 
COBE, large-scale peculiar velocities, and the power spectrum inferred from redshift surveys of 
IRAS galaxies. 



II. PARTICLE PHYSICS MODELS 



There are a number of ways in which massive pseudo-Nambu-Goldstone bosons with the 
requisite mass scales discussed above may play a role in particle physics models. In this section, 
we schematically outline only a few of them. The basic idea is to build a model with a global 
symmetry spontaneously broken at a large mass scale / ~ Mpi, which gives rise to one or more 
massless Nambu-Goldstone bosons. There are then several ways to introduce explicit breaking of 
(some or all of) the global symmetry at the scale A ~ Mqut, resulting in potentials for the would- 
be Goldstone modes. Ideally, the lower scale emerges dynamically, so that no small parameters 
are introduced. 

The most familiar example of a pseudo-Nambu-Goldstone boson in nature is the pion. Here, 
the global chiral symmetry is spontaneously broken by quark condensates at the QCD scale, 
(qq) ~ Aq CD ~ (100 MeV) 3 , and explicitly broken by quark masses, m u ~ m d ~ 10 MeV. In 
the case of the pion, these two scales are close together (they differ by a factor of about ten), 
so the pion gains a mass comparable to the QCD scale, m 2 ~ m q (qq) / f 2 ~ (100 MeV) 2 . By 
contrast, in invisible axion models [7], the scales of spontaneous and of explicit symmetry breaking 
are separated by many orders of magnitude: the spontaneous symmetry breaking scale /pq is 
elevated close to the GUT scale, while the explicit breaking scale is ~ Aqcd- The resulting 
hierarchy of scales yields a very light axion, m 2 ~ rn q (qq) / fpQ~, for example, m a ~ 10~ 5 eV for 
fpQ — 10 12 GeV. For the PNGB inflaton, we will be interested in models with a relatively modest 
hierarchy between the spontaneous and explicit global symmetry breaking scales, A// ~ 10 -4 . 
Such a ratio of scales is intermediate between the case of the pion (A// ~ 0.1) and the invisible 
QCD axion (A// - 10" 13 ). 

A. PNGBs from Condensates 

In this section, we illustrate how such an intermediate mass hierarchy can arise. We consider 
an action that contains coupled scalar and fermion fields and exhibits a chiral U(l) symmetry. 
Spontaneous breaking of the chiral symmetry takes place at energy scale / (for inflation, / ~ M p i); 
massless Nambu-Goldstone bosons arise at this scale. We illustrate an additional feature that may 
be attractive although not necessary to our model: if the scalar field couples non-minimally to 
gravity, it may dynamically generate Newton's constant at this scale (induced gravity) [12]. Next, 
we discuss several ways in which the symmetry can be explicitly broken at a lower energy scale 
~ A (for inflation, A ~ 10 _4 M p z). At this scale, the Nambu-Goldstone boson acquires a mass, 
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in a manner similar to the axion or schizon [13] (although at higher mass scale). We focus on 
axion-like scenarios, in which a gauge group becomes strong at the scale ~ A. We briefly discuss 
how this may arise in technicolor models and then, in somewhat more detail in Sec. IIB, in 
superstring models. 



1) Spontaneous Symmetry Breaking 

Taking our cue from the axion [14], we first describe a simple model which implements the 
mechanism described above. Consider the fundamental action for a complex scalar field <I> and 
fermion ip, coupled to gravity: 

S = j d 4 x^g~ [g^d^d,,® - V($*<S>) - £<f>*$R + i^d^ip - (h$ L i> R <f> + h.c.)] (2.1) 

where R is the Ricci scalar, and ip(R,L) are respectively right- and left-handed projections of 
the fermion field, 4>(r,l) = (1± 7 5 )V'/2- This action is invariant under the global chiral U(l) 
symmetry: 

ipL e iaj2 ^L , *I>r e~ ia /^ R , $ e ia $ , (2.2) 
analogous to the Peccei-Quinn symmetry in axion models. 

We assume the global symmetry is spontaneously broken at the energy scale / in the usual 
way, e.g., via a potential of the form 

F(|$|) = A^*$-^ , (2.3) 

where the scalar self-coupling A can be of order unity. The resulting scalar field vacuum expecta- 
tion value (vev) is (<J>) = je^l f / y/2. 

In this model, spontaneous symmetry breaking dynamically generates Newton's constant for 
Einstein gravity [12]. At scales below /, the non-minimal coupling of the scalar field to the 
curvature induces the canonical Einstein Lagrangian, £(<!>*<I>)i? = (^f 2 /2)R = R/16ttG, if the 
coupling £ satisfies 

(24) 

Since inflation requires / ~ Mpi, the above relation holds for £ of order unity, a natural value 
for this dimensionless coupling. We note that generation of the Planck scale in this way is not a 
necessary ingredient of the models discussed below: since inflation takes place in the <p direction, 
after reaches its vev, we could simply replace the non-minimal coupling term in Eq.(2.1) with 
the usual Einstein Lagrangian. On the other hand, since the mass scale / must be comparable 
to Mpi for successful inflation, it is natural and economical to tie it directly to the gravitational 
scale. Since the gravitational sector is canonical once the temperature of the universe drops below 
the scale /, we assume ordinary Einstein gravity from now on. 

Below the scale /, we can neglect the superheavy radial mode of $ (m ra di a i = A 1 / 2 / ~ Mpi) 
since it is so massive that it is frozen out. The remaining light degree of freedom is the angular 
variable (ft, the Goldstone boson of the spontaneously broken U(l) (one can think of this as the 
angle around the bottom of the Mexican hat described by eqn. (2.3)). We thus study the effective 
chiral Lagrangian for (f>: 

£eff = \d^4> + i^d^iP - (m ^Re i<t>/f + h.c.) . (2.5) 
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Here the induced fermion mass mo = hf/-\/2; for example, for values of the Yukawa coupling 
10 -3 < h < 1, the fermion mass is in the range Mqut < m o < M p [. The global symmetry is now 
realized in the Goldstone mode: C e ff is invariant under 

V>l - e ia ^ L , ift R - e" iQ / 2 Vfl , (ft - <ft + af . (2.6) 
At this stage, (ft is massless because we have not yet explicitly broken the chiral symmetry. 

2) Explicit Symmetry Breaking 

Several options exist for explicitly breaking the global symmetry and generating a PNGB 
potential at a mass scale ~ A several orders of magnitude below the spontaneous symmetry 
breaking scale /. In a class of .^-symmetric models studied by Hill and Ross [13], one adds a bare 
fermion mass term mi'ipL'^R to £ e ff, which presumably arises from another sector of the theory 
(just as quark masses in QCD are generated in the electroweak sector). The combination of terms 
involving mo and mi generates a 1-loop potential for (ft of the form (1.1), with A 2 ~ momi; a 
synopsis of these 'schizon' models is given in refs. [13,15]. 

For the rest of this discussion, we focus on the simplest mechanism for explicit symmetry 
breaking, by analogy with the QCD axion: dynamical chiral symmetry breaking through strongly 
coupled gauge fields. Suppose the gauge symmetry of the effective theory below the scale / ~ Mpi 
is a product group, G\ x G2, where G\ is a standard grand unified group {e.g., E 6 or SU(5)) 
which spontaneously breaks down to the standard model at some scale Mqut- In other words, 
G\ describes the physics of ordinary quarks and leptons (and their heavier brethren) while G2 
might describe a 'hidden sector'. At the G\ unification scale, the G\ gauge coupling is small 
(perturbative unification). On the other hand, let G2 be an asymptotically free non-abelian 
gauge theory which becomes strongly interacting at a scale k comparable to the GUT scale. In 
addition, we assume that tft transforms non-trivially under G2 (ift carries GVcolor'). Starting 
with a perturbative G2 gauge coupling at the Planck scale, 02 (Mpi) = g\{Mpi)jAm (which is, 
say, comparable to ai(Mpi)), the scale k emerges from the renormalization group, 

where the renormalization group constant 60 determines the lowest order term in the expansion of 
the /3-function of G2, ^2(9) = —bog^/i^) 2 — .... For example, for G2 = SU(N) and no light matter 
fields with G2 charge, then 6 = 3A; if there are N matter fields (one generation) with masses 
m < k in the fundamental representation of G2, then 6 = 2 A. For reasonably large groups, 
and therefore large 60, the gauge coupling can run sufficiently fast to generate k ~ Mqut- As 
examples, for a 2 (-^p«) = 1/30 and G2 = SU(5) we find k ~ 3 x 10 14 GeV if there are no light 
(m < Mpi) fermions transforming under G2; on the other hand, with N light fermions, the same 
value of k arises for the larger group G2 = SU (9). 

Since ift is charged under G2, we expect chiral dynamics to induce a fermion condensate, 
(■tfttp) ~ ft 3 . (We assume the condensate can be rotated to be real; the extra phase it involves 
is irrelevant for our discussion). From eqn.(2.5), the condensate explicitly breaks the global 
symmetry, giving rise to a potential for the angular PNGB field (ft, 

V{(ft) = Re[m (^Vii)e^ // ] = ™o« 3 cos(</>//) . (2.8) 

This has the form of eqn.(l.l), with A 4 = m K 3 = /i/k 3 /\/2. For inflation, we require A ~ Mqut- 
Such an energy scale can arise in at least two ways: (i) mo ~ k ~ Mqut', this requires the 
Yukawa coupling h ~ 10~ 4 , or (ii) m ~ Mpi, h = O(l), and k is slightly below the GUT scale, 
k ~ 10 _1 Mg;7t- We indicated above that the running of the coupling constant for group G2 
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may indeed provide such a value for k. For this second choice of parameters, we do not need to 
introduce any small coupling constants in the fundamental Lagrangian near the Planck scale: the 
small ratio A// emerges dynamically and is "strongly natural". 

Although this model may be cosmologically appealing, we do not want to propose a new 
strongly interacting gauge sector in particle physics solely to generate an inflaton potential. Hap- 
pily, there is well-founded particle physics motivation for an additional gauge group which becomes 
strong at the GUT scale, and this idea has a distinguished history in the particle physics literature. 

One possibility is that G2 is a technicolor group, and that ip carries both G\ and G2 charge. 
[In this case, (j> can couple through ^ — — ^-triangle diagram to ordinary particles (e.g., 
gluons and photons); this may be advantageous in that it leads to reheating of the 'ordinary' 
sector of quarks and leptons. We thank S. Dimopoulos for making this point to us.] Then, one 
must introduce a source for spontaneous breaking of the standard G\ GUT group. Here, one 
may contemplate two possibilities: If the (iftip) condensate is a G\ singlet (this may happen even 
though ip carries G\ charge), then G\ must be broken by the usual Higgs mechanism or some 
equivalent. Alternatively, if the condensate (i)ip) is G\ -non-singlet, it can spontaneously break 
G\ at a scale k ~ Mgut, by analogy with technicolor models. If it can be implemented, the latter 
choice would be most economical: a single mechanism would give rise to both GUT symmetry 
breaking and inflation, and the only fundamental scalar (<&) in the theory has Planck mass. This 
value of the mass for a scalar is natural, and in principle no small parameters would need to be 
introduced in the theory. 

B. A Superstring Model: "Supernatural" Inflation 

A second motivation for a gauge group which becomes strongly interacting at the GUT 
scale comes from superstring theory [16]. In these models, the gauge symmetry of the effective 
supergravity theory below the Planck scale is again a product group, G\ x G2, where G\ = Gqut 
contains the standard model and G2 describes the hidden sector; for example, in the original 
heterotic string model, G\ x G2 = Eg x E' 8 . 

In the effective field theory arising from superstrings, an important role is played by the 
complex scalar field S. The real part of this field, ReS", is the dilaton; the imaginary part ImS* 
is the 'model-independent axion'. In string theory, the value of the dilaton determines the string 
coupling constant g s through the relation [17] (Re(S)) = l/g^- Since g s is related by factors of 
O(l) to the gauge couplings g a (M p i) of the effective field theory at the Planck scale (a labels the 
gauge group G a ), the dilaton expectation value determines gauge couplings as well. In particular, 
a value for the dilaton in the range (ReS) = ReSo ~ 1.5 — 2.5 yields a phenomenologically viable 
G\ gauge coupling at the GUT scale, q.\{Mgut) ~ 1/45. If this theory is to have predictive 
power and time-independent constants of nature, one would expect that the dilaton potential 
U(Re£) has a minimum in this range. In perturbation theory, the dilaton (and axion) potential 
V(S) is protected by supersymmetry: if supersymmetry is unbroken at tree level, then V(S) 
vanishes to all finite orders in perturbation theory, leaving the gauge couplings indeterminate 
[18]. However, if the hidden sector G2 is an asymptotically free non-abelian group, it will become 
strongly interacting, leading to condensation of gauginos (fermion supersymmetric partners of the 
gauge bosons) at a scale (AA) ~ M^e -8 "" 2 / b °9l( M Pi) _ Through the relation between ReS" and the 
string coupling constant above, this corresponds to a non-perturbative potential for the dilaton 
of the form V(S) oc e~ cS . As a consequence, the imaginary part of the field, the axion partner 
of the dilaton, obtains a potential of the form (1.1), V(ImS) oc cos(ImS'). Our interest in this 
scenario derives from the fact that the 'model-independent' axion could in principle play the role 
of the inflaton in natural inflation [19]. 

Although nonperturbative effects in the hidden sector can generate a potential for the dilaton, 
an exponentially falling potential clearly will not by itself stabilize the dilaton in the desired range 
noted above: instead ReS runs away to infinity, yielding a free string theory. Additional physics is 
needed to help pin the dilaton at the appropriate minimum; we describe this further below. Here, 



8 



we mention that a second important role in particle physics of hidden sector gaugino condensation 
is that it might break supersymmetry (SUSY) [9]. If the condensate breaks supersymmetry in 
the hidden sector at the scale (AA) ~ Mqut ~ 10 14 GeV, then SUSY is broken in the observable 
sector at the scale Msusy ~ Mq UT /M 2 [ ~ TeV. SUSY breaking at this scale would protect 

the small Higgs mass and alleviate the heirarchy problem. Thus, the factor e -1 ^ 2 in the scale 
of gaugino condensation might lead to a large gauge hierarchy. (It is also possible that SUSY 
breaking arises from some other mechanism.) 

The first attempts to implement these ideas in the G\ x G2 = E 8 x E' 8 heterotic string 
theory relied on the hidden E' s sector becoming strongly interacting, and generating gaugino 
condensation, at a scale comparable to the GUT scale [9]. As noted above, in this case the 
gaugino condensation-generated potential for KeS decays exponentially for large values of the 
dilaton field. Attempts were made [9] to stabilize the dilaton by combining gaugino condensation 
with a term arising from the expectation value of the antisymmetric tensor field {H^x). However, 
quantization conditions on the vacuum expectation value of this field [20] require it to be of order 
unity in Planck units, implying that the resulting potential V(S) only has a minimum where KeS 
is small (well below the desired range above), i.e., where g s is large. As a result, the string theory 
would be strongly coupled, and the whole framework of perturbative calculations in the effective 
field theory would be unreliable [21]. 

Recently, this problem has been reconsidered by Krasnikov [22], Casas, et al. [23], and 
Kaplunovsky, Dixon, Louis, and Peskin (hereafter KDLP) [24], in the context of string models 
where the hidden sector G 2 is itself a product of two or more gauge groups. They found that 
the combined effect of gaugino condensates in multiple hidden groups can generate a dilaton 
potential with a weak-coupling (small perturbative g s ) minimum. In some cases, supersymmetry 
also appears to be broken at the requisite scale (L. Dixon, private communication; Kaplunovsky, 
unpublished). Here we will briefly study the axion potential generated in these multiple gaugino 
condensate models and explore its suitability for inflation. 

The effective Lagrangian for the dilaton S can be written 

M 2 

Cs = Msf¥f d " sd ^ ~ m sl (2 ' 9) 

where V is the effective potential generated by gaugino condensation. In string theory, the Planck 
scale is derived from the fundamental string tension a' via M Pl = IQtt /g 2 s a'. At tree level, the 
gauge coupling of group G a is then g a = g s j y/k^, where k a is the level of the Lie algebra of G a 
(a small integer). Thus, at tree level, we have (ReS) = ReS'o = [4irk a aGUT(Mpi)]~ 1 ; assuming 
G\ (Gqut), which contains the standard model, is at level one {k\ = 1), the phenomenologically 
acceptable value of the GUT gauge coupling, aGUT(Mpi) ~ 1/20 — 1/30, requires that the dilaton 
VEV be in the range ReS'o = 1-5 — 2.5, as noted above. As KDLP show, this large a value of the 
dilaton expectation value can be obtained with a hidden group structure Gi = SU{N\) x Sf^A^), 
provided that the expression [(&i/Aq) — (A^/A^)] -1 is large; e.g., for their 'best case', k\ = k 2 = 1 
and Aq = 9, N 2 = 10 (see below). In what follows, for simplicity, we shall follow KDLP in taking 
G2 to be a product of two SU (N) groups. 

Following KDLP and ignoring gravitational and subleading 1/N corrections (that is, con- 
sidering a global SUSY model with large hidden gauge groups), the effective dilaton potential 
is 

V (S,S*) = ^(S + S^l J>(AA> a | 2 (2.10) 

PI a 

where subscript a = 1, 2 now refers to the hidden gauge group G2 a . When the coupling constant 
of group G2 a = SU(N a ) becomes strong, the resulting gaugino condensate is 

(AA) a = N a vM? en e ie » exp 



24x z k a S + §A a 



'J0,a 



(2.11) 
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Here, the renormalization mass scale (at which the effective Lagrangian is defined) is taken to be 

M 2 ren = £ ; p 2 = ^= = (0.216) 2 , (2.12) 

where 7 is Euler's constant and a' = 8w(S + £*)/ME is the inverse string tension; v is an N- 
independent constant of order unity; 9 a = 2irm/N a , with m integer, is an arbitrary discrete phase 
reflecting the iV a -fold degeneracy of the vacuum states of the theory; 60, a is the renormalization 
group constant for group G 2a ; and A a is the threshold renormalization factor of order N a . A a , 
which in general can be a function of the moduli fields T l , enters the coupling constant to one- loop 
order via 

: l US + S*)-^-Ao g ^ + ^. (2.13) 



gl(n) 2 ay ' 16vr 2 & ^ 16vr 2 

Decomposing the dilaton into real and imaginary parts, and assuming no charged fermions in the 
hidden groups (i.e., taking & ,a = 3JV a ), we have 



v 2 M 4 

ns) = 5xl0 -»^ 



( 16^ 2 fc 1 R.eS + A 1 \ / 167r 2 Jc 2 ReS + A 2 \ 



2k 1 k 2 N 1 N 2 e(-^(^ + ^ S -^ + ^)) cos ( 8 vr 2 (jl- - ^) ImS + Sej 



(2.14) 



where 80 = 9 2 — 9 1 . 

To study inflation, it is preferable to work with scalar fields that have canonical kinetic terms 
in the Lagrangian. From eqn.(2.9), the kinetic Lagrangian for the real and imaginary components 
of the dilaton is not of this form, since 

M 2 

Chin = ™ 2 (dpRieSd»BeS + ^ImSc^ImS) (2.15) 



Thus the canonically normalized real component is taken to be [19] (f>R = —Mpi ln(ReS')/v / 167r. 
In general, the real and imaginary parts of S are interdependent, and one should follow the 
coupled evolution in the two-dimensional field space. For simplicity, to focus on the imaginary 
component, the model-independent axion, we shall assume the real component reaches its VEV, 
(ReS) = ReSo, well before the imaginary part does; in a chaotic scenario in which the field is 
initially randomly distributed, this will always be true in some regions of space. (Note that, near 
the Planck time, S will drop out of thermal equilibrium and its potential will be dynamically 
negligible; under these conditions, we expect no special initial value for S to be preferred.) In 
that case, we can define the canonical axion field, 

(f> a = MpilmS/VmrReSo (2.16). 

However, from eqns. (1.1) and (2.14), we have (f> a /f = 87r 2 ImS'[(£;i/iVi) — (A^/iNfo)]. Combining 
these two expressions, we find the equivalent global symmetry breaking scale 

/= ^ (h_h_y 1 

As we will see in Sections III and IV, the phenomenologically acceptable range for / is / > 0.3Mpj. 
From eqn.(2.14), the potential for the real part of the dilaton is minimized at 



P o Wfci k 2 
ReS ° ~ 8^ I m ~ W 2 



ko \No 2 



(2.18) 
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where 



A2 _ Ai 

N 2 iVi 



(2.19) 



is the difference in threshold renormalization factors. We thus find 



Mpi r (k x 




1 -1 



/ = 



(2.20) 



In order to achieve an acceptably large value for ReSo, KDLP choose, e.g., N± = 9, N 2 = 10, with 
k\ = k 2 = 1; larger values of iV a are excluded because the size of the hidden sector is constrained 
by the total Virasoro central charge available. With this choice, to obtain ReSo > 1.5 requires 
5 A > 2.8, which implies f/M P i =< 0.11 ~ 1/V24^. If this upper limit is saturated, a sufficiently 
long epoch of slow-rollover inflation can occur (see section III), but the reheat temperature is 
unacceptably low and the density perturbation spectrum has too much power on large scales. On 
the other hand, for 8 a = 1, we would have f /Mpi ~ 0.36, which yields a viable inflationary model 
with an interesting fluctuation spectrum. However, as eqn.(2.18) shows, this value of <5a would 
require larger groups, e.g., N 1: N 2 = 16, 17, to achieve ReSo > 1-5, and this violates the central 
charge limit (however, see comment below). 

For these models, we can read off the effective scale A, defined in eqn.(l.l), from eqns.(2.14) 
and (2.18); A is determined up to a constant of order unity (the factor v) by the values of k a , N a , 
and A„. For the SU(9) x SU(W) example, taking 5 A = 2.8, A1/JV1 = -A 2 /N 2 = -1.4, which 
corresponds to ReS = 1.5 and / = 0.11M P ; ~ Mpij^fl^K, we find A = 6 x 10 _ V/ 2 M PJ = 
8 x lO 14 !; 1 / 2 GeV, in the right vicinity for generating an acceptable density fluctuation amplitude 
(even though, as noted above, this value of / leads to an unacceptable fluctuation spectrum-see 
section IV). This is a pleasing feature of these models: the same physics which sets the condensate 
scale to be of order Mqut (~ Mp;exp[— 8ir 2 /g 2 bo]) fixes A to approximately the same scale. 

From our perspective, the interesting result here is that a string model designed to yield a 
phenomenologically plausible particle physics scenario, in particular a large gauge hierarchy and 
possibly supersymmetry breaking near the weak scale, implies values for the PNGB parameters 
/ and A for the model-independent axion which are quite close to those needed for successful 
inflation. Furthermore, as suggested in [23,24], with the inclusion of charged matter fields in one 
of the hidden groups, it is possible that the value of 5 A required to fix ReSo could be reduced 
from ~ 3 to ~ 1, generating a sufficiently large value of / for inflation. 

We end this subsection with several caveats about the treatment given here. First, as men- 
tioned above, we have reduced a two-dimensional problem to a one-dimensional one by assuming 
the dilaton is already pegged to its expectation value during the evolution of the axion. Although 
this will be accurate for some region of parameter space, and for chaotic initial conditions in 
some regions of the universe, in general one should treat the full two-dimensional problem. In 
particular, the possibility of inflation in the dilaton direction deserves study (the potential for the 
canonical dilaton contains terms of the form exp[— ae~ b< ^ R \). Second, due to 2-loop running of 
the gauge coupling, the prefactor of the exponential in eqn.(2.11) actually contains an additional 
factor of S [24]. This gives rise to an overall multiplicative factor of SS* on the right side of 
eqn. (2.14), modifying the dependence of the potential on the axion field from a pure cosine. 
This could have interesting consequences for the cosmological evolution of the model-independent 
axion. Third, in this discussion we have assumed that the dominant non-perturbative effects in 
string theory arise at the level of the effective supergravity Lagrangian. It has been suggested [25] 
that some inherently stringy non-perturbative effects at the Planck scale are only suppressed by 
a factor exp(— 2n/g) as opposed to the field theory factor exp(— 8TT 2 /g 2 ). If such stringy effects 
contribute to the effective dilaton potential, they could substantially modify this effective field 
theory analysis. 

C. Alternatives 
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In the preceding subsections, we have outlined two particle physics models which incorporate 
a PNGB with the requisite parameters for inflation. Clearly there are further possibilities [13,15]. 
For example, one can imagine doing away with fundamental scalars altogether, and having the 
PNGB arise as an effective field. One choice would be a composite PNGB built from a fermion 
condensate, in analogy with composite axion models [14] and the pion. A second possibility, 
recently discussed by Ovrut and Thomas [26] , builds on the existence of instantons in the theory 
of an antisymmetric tensor field B^ v (recall that such a tensor field arises, e.g., in superstring 
theory.) Defining the field strength 

H n»x = Qngn* + g*B»" + d v B x » , (2.21) 

the action for this theory is 

S = ^2 / d A xH^ x H^ x , (2.22) 

with the resulting equation of motion d^H^ uX = 0. As Ovrut and Thomas note, the theory 
(2.22) has pointlike, singular instanton solutions, analogous to Dirac monopoles in electromag- 
netism; evaluating their contribution to the partition function, the resulting effective action can 
be expressed in terms of an effective mean scalar field 4> as 



d 4 x 



\e%m ~ 2Ve"^ 4 ) 1/3 («V) 1/3 f (l ± cos (|)) 



(2.23) 



where a is a number, and / is a mass scale characterizing the instanton solutions. Clearly this 
is of the form (1.1) and, for / ~ Mpi, (2.23) is another potential candidate model for natural 
inflation. In a variant of these models, the tensor field can be coupled to a fundamental real 
scalar field u with the symmetry-breaking potential V(u) = (A/4!)(it 2 — 6ra 2 /A) 2 . This also leads 
to a potential of the form (1.1) for the associated scalar mean field theory; for / ~ m ~ Mpi 
and A ~ 10 -4 , one finds [26] A ~ 10 16 GeV, as desired for successful inflation. In both of these 
models, as in the string model of the previous subsection, the effective scale A is small compared 
to / due to the exponential (instanton) suppression factor. This is the origin of the hierarchy 
required for the generation of acceptably small density fluctuations in inflation. The advantage 
of these models is that this hierarchy does not need to be put in by hand. 



D. Other Issues 

Before leaving this survey of model-building, we note recent work drawing attention to the 
fact that global symmetries may be explicitly broken by quantum gravity effects [27,28] (e.g., 
wormholes and black holes). If such effects are characterized by the Planck scale, they may 
induce non-renormalizable higher-dimension terms in the low-energy effective Lagrangian for <I> 
(see eqn.2.3), of the form 

Veff(*) = 9mn ^2m+n-4 " ( 2 - 24 ) 

The coefficients g mn introduced here should not be confused with the gauge and string couplings 
discussed above. Terms with n / explicitly break the global U(l) symmetry of eqn.(2.2). Taking 
9mn = |5'mn|exp(i(5 mn ), the induced PNGB potential is a sum of terms of the form 

(j, \ 2m+n / , \ 

M Pl co S rj + 5 mn \ . (2.25) 

Therefore, for n / 0, the effective explicit symmetry breaking scale is 

1m-\-n 

K ef f = \g mn \ 1 ^M Pl U [ -\ 4 , (2.26) 
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where non-renormalizable terms have dimension 2m + n > 5. Since PNGB inflation requires 
/ > 0.3Mpi and A ~ Mqut, the coefficients g mn of these terms must be relatively small; for 
example, for the dimension 5 term, g < 10 -14 is required. (The upper limit on g mn is relaxed for 
higher dimension terms.) 

Naively, this effect appears to lead us back to the same difficulties this inflation model was 
meant to solve, namely, a small dimensionless constant of order 10~ 14 appearing in the Lagrangian. 
However, it is worth making several remarks about this problem. First, a caveat: in the discussion 
above (and in refs.(28)), it was implicitly assumed that the coefficients g mn are 'naturally' of order 
unity. However, in the absence of a solvable quantum theory of gravity, these coefficients cannot be 
reliably calculated. In model wormhole calculations, one must introduce a cutoff scale [i <C Mpi, 
in which case such effective operators are proportional to the tunneling factor ~ exp(— M Pl /fi 2 ). 
Thus, in the regime in which one can calculate, the coefficients g mn are highly suppressed; the 
assumption that they are not at all suppressed depends on an uncertain extrapolation of the 
cutoff scale to the Planck scale. In addition, there may be other effects which enter to suppress 
these terms. In particular, in the axion model studied in more depth in ref.[27], wormhole effects 
are effectively cut off at the symmetry breaking scale /, leading to an exponential suppression 
~ exp(— Mpi/f) in the wormhole induced axion potential. Second, even supposing such terms 
are in principle unsuppressed (all g mn of order unity), there are ways in which they could be 
evaded. For example, for a large gauge group (as contemplated above), a global symmetry may 
automatically be present, due to the gauge symmetry and field content of the theory, preventing 
terms up to some relatively large value of 2m + n; in the present case, this would require that 
all terms up to 2m + n ~ 25 be forbidden. Alternatively, if the field (f> is an effective field which 
arises below the Planck scale, as suggested above, such explicit symmetry breaking terms can 
be forbidden by a local symmetry of the underlying theory, as in the superstring example of 
the previous subsection. Alternatively, as in the antisymmetric tensor model of section C above, 
the (j) field may be unrelated to a global symmetry. Therefore, while the arguments of [28] are 
provocative, there are many examples of particle physics models which give rise to potentials of 
the form (1.1), with the requisite mass scales for inflation, which evade them. 



III. COSMIC EVOLUTION OF THE INFLATON FIELD 



With these models as theoretical inspiration, we turn now to the cosmological dynamics 
of an effective scalar field theory with a potential of the form (1.1) below the scale /. For 
example, in the model of Sec. II. A, / is the global spontaneous symmetry breaking scale, and 
4> describes the phase degree of freedom around the bottom of the Mexican hat potential (2.3); 
in other models, however, the picture may differ. To successfully solve the cosmological puzzles 
of the standard cosmology, an inflationary model must satisfy a variety of constraints, including 
sufficient inflation (greater than 60 e-folds of accelerated expansion) for a reasonable range of 
initial conditions; sufficiently high reheat temperature to generate a baryon asymmetry after 
inflation; and an acceptable amplitude and spectrum of density fluctuations. In this section we 
explore these constraints analytically and numerically for potentials of the form (1.1). 

The 4> interaction cross-sections with other fields are generally of order a ~ l// 2 , so its 
interaction rate is of order T 3 // 2 . Comparing this with the expansion rate H ~ T 2 /Mpi, 

we see that the scalar inflaton field thermally decouples at a temperature T ~ f 2 /Mpi ~ /. 
We therefore assume 4> is initially laid down at random between and 2irf in different causally 
connected regions. (This is the simplest but by no means only possible initial condition.) Within 
each Hubble volume, (i.e., ignoring spatial gradients-see below) the evolution of the field is then 
described by the classical equation of motion for a homogeneous field 4>(t), 

4> + SH<j) + r<j> + V'(<f>) = , (3.1) 
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where T is the decay width of the inflaton, and the expansion rate H = a/ a is determined by the 
Einstein equation, 



H 1 



8ir 



v(4>) + -4> 2 



For completeness, it is also useful to have the second order Friedmann equation, 



8tt 
3Mf t 



<j> 2 - V(<f>) 



(3.2) 



(3.3) 



In Eqns. (3.2-3), we have assumed that the scalar field dominates the stress energy of the universe; 
this will hold starting near the onset of inflation. 

In the temperature range A < T < f, the potential V(4>) is dynamically irrelevant, because 
the forcing term V'{4>) in Eq.(3.1) is negligible compared to the Hubble-damping term. (In 
addition, for axion-like models in which V(4>) is generated by non-perturbative gauge effects, 
A — > as T/A — > oo due to the high-temperature suppression of instantons [14].) Thus, in 
this temperature range, aside from the smoothing of spatial gradients in (f> (see below), the field 
does not evolve. Finally, for T < A, in regions of the universe with <p initially near the top 
of the potential, the field starts to roll slowly down the hill toward the minimum. In those 
regions, the energy density of the universe is quickly dominated by the vacuum contribution 
(V(4>) ~ 2A 4 > p ra( i ~ T 4 ), and the universe expands exponentially. Since the initial conditions 
for (f> are random, our model is closest in spirit to the chaotic inflationary scenario [29]. In 
succeeding subsections, we study this evolution in more detail. 



A. Standard Slow-Rollover Analysis 

In this subsection, we recapitulate the analytic treatment of PNGB inflation given in Paper 
I. A sufficient, but not necessary, condition for inflation is that the field be slowly rolling (SR) 
in its potential. Therefore, by analyzing the conditions for, and number of e-foldings of, inflation 
in the SR regime, we should be at worst underestimating the true number of inflation e-folds. 
The field is said to be slowly rolling when its motion is overdamped, i.e., (f> <C 3H(p, so that the 
(f) term can be dropped in Eq. (3.1) (n.b., we assume r <C H during this phase). It is easy to 
show that in general this SR condition is a sufficient condition for inflation. First, from the scalar 
equation of motion (3.1), the defining SR condition implies that <j) 2 <C 2V(4>). On the other hand, 
the universe is inflating if the Robertson- Walker scale factor a(t) is accelerating, a > 0; from Eqn. 
(3.3), this requires 4> 2 < V. Thus, if the SR condition is well satisfied, we are guaranteed to be 
in an inflationary epoch. The converse is not necessarily true: inflation can occur even when the 
field is not slowly rolling. However, we will see in subsequent sections that, for this potential, if / 
is larger than about Mp\ / \/24it, the SR epoch is roughly coincident with the inflationary epoch. 

Hereon, for the purposes of numerical estimates, we shall assume inflation begins at a field 
value < (f>i/f < 7r; since the potential is symmetric about its minimum, we could just as easily 
consider the case ir < (f>i/f < 2tt. For the potential (1.1), the SR condition implies that two 
conditions are satisfied: 



and 

V'{<t>)M Pl 



v{4>) 



l + cos(0//) Mpi 



From Eqns. (3.4), the existence of a broad SR regime requires / > Mp\j \/487r (required below 
for other reasons). The SR epoch ends when (f> reaches a value fa, at which one of the inequalities 
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(3.4) is violated. In Fig.l, we show fa/ f as a function of f /Mpf, as / grows, fa/ f approaches 
the potential minimum at it. For example, for / = Mpi, (j) 2 /f = 2.98, while for / = Mpi/ 
fa/ f = 1-9. For / > 0.3Mpi which, as we shall see below, is the mass range of greatest interest, 
the two inequalities (3.4a) and (3.4b) give very similar estimates for fa. For simplicity, we can 
then use (3.4b) to obtain 

^2arctan^^|^j (/ > 0.3M PI ) . (3.5) 

Once <p grows beyond fa, the field evolution is more appropriately described in terms of oscillations 
about the potential minimum, and reheating takes place, as described below. We note that the 
expansion of the universe (the 3H(p term in Eq. 3.1) acts as a strong enough source of friction 
that the field is not able to roll through the minimum at irf and back up the other side sufficiently 
far to have any further inflationary period. 

To solve the standard cosmological puzzles, we demand that the scale factor of the universe 
inflates by at least 60 e-foldings during the SR regime, 



f t2 Sir f 4 " 2 V(6) 

N e{ fa,fa, f) , H a 2 /a^^ H*=j^] ^ 



16tt/ 
M Pl 2 



2 



In 



Mfam 

sm(fa/2f)\ 



> 60 . (3.6) 



Using Eqns. (3.4,5) to determine fa as a function of /, the constraint (3.6) determines the 
maximum initial value {fa^ ax ) of fa consistent with sufficient inflation, N e ((j)™ ax , fa, f) = 60. For 
/>0.3M P; , 



M 



2 



sin(<^r72/)^ 1 + ' exp(-^^) . (3.7) 



1/2 / 15M, 



2 



1 pi 



48tt/ 2 J 1 \ 4tt/ : 

The fraction of the universe with fa 6 [0, ^>™ ax ] will inflate sufficiently If we assume that fa 
is randomly distributed between and nf from one Hubble volume to another, the a priori 
probability of being in such a region is P = fa^ ax /irf '. For example, for / = 3Mpi, Mpi, M p i/2, 
and Mpij^p2AK, the probability P = 0.7, 0.2, 3 x 10 -3 , and 3 x 10 -41 . The initial fraction of 
the universe that inflates sufficiently drops precipitously with decreasing /, but is large for / 
near Mpi. This is shown in Fig. 1, which displays \og{(f>™ ax / /) = 0.5 + logP and fa/ f. These 
considerations show that for values of / sufficiently near Mp\, sufficient inflation takes place for 
a broad range of initial values of the field <fi. We note that these constraints do not determine the 
second mass scale A. 

According to some inflationists, the discussion above of the probability of sufficient inflation 
is overly conservative, since it did not take into account the extra relative growth of the regions 
of the universe that inflate. After inflation, those initial Hubble volumes of the universe that did 
inflate end up occupying a much larger volume than those that did not. Hence, below we will 
also compute the a posteriori probability of inflation, that is, the fraction of the final volume of 
the universe that inflated. 



B. Numerical Evolution of the Scalar Field 

In this section, we expand upon the results of the preceding subsection by numerically inte- 
grating the equations of motion. This yields a more accurate estimate of the time (or field value) 
when inflation ends and the amount of inflation that takes place, as a function of the mass scale 
/ and the initial value of the field fa . 
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First we rewrite Eq. (3.1) in terms of more useful variables. As a dimensionless time variable, 
we use the number of e-foldings of the scale factor, 



dn = Hdt, (3.8) 

and hence d/dn = H~ 1 d/dt. We also define the dimensionless field value, field 'velocity', and 
mass ratio 

y = faf ; v = dy/dn ; 7 = 3M^/8vr/ 2 . (3.9) 
Then we can write Eqs. (3.1) and (3.2) as 



dv 
dn 



3/2 



7 tan(y/2) - 3v 1 - v 2 /2-f - ujv 1 - v 2 /2^ "' "7 1 / 2 [l + cos y] 1/2 , (3.10) 



where oj = Tf/A 2 contains the effects of dissipation. For the purpose of numerically calculating 
the evolution of the field, we will assume as in the previous section that this dissipation term is 
negligible. In this approximation, Eq. (3.10) depends only on the shape of the potential and on 
7, i.e., on the ratio f/Mpi, and not explicitly on A. 

In order to solve Eq. (3.10), we must specify two initial conditions: the initial values of the 
field and its time derivative. We allow the initial field value y\ = fa/ f to range over the interval 
to 7r, and take the initial velocity to be v\ = fa/Hf = 0. The assumption of zero initial velocity 
is the one usually made in discussions of inflationary models. However, in the course of smoothing 
out gradients or due to randomness in the initial conditions, we expect the field to acquire an 
initial 'Kibble' [30] velocity at the temperature T ~ A such that its kinetic energy is comparable 
to the potential energy ~ A 4 . Naively, this velocity effect could delay or even prevent the onset of 
inflation. This problem has been studied previously in the context of new and chaotic inflationary 
models and has been shown to be potentially problematic for new inflation [31]. Initial velocities 
in the context of the present model have been studied numerically by Knox and Olinto [32]. They 
find that, due to the periodic nature of the potential, the effect of initial velocities is merely to 
shift, but not change the size of, the phase space of initial field values which lead to at least 
60 e-folds of inflation. That is, as fa is increased from zero, the value of fa at which inflation 
begins is shifted, but the fraction of initial field space which inflates is approximately invariant. 
Therefore, for models of the form (1.1), we lose no generality by assuming fa = 0. Given these 
initial conditions, we solve the equation of motion (3.10) numerically. The resulting solution 
y(n,yi) provides the value of the 0-field after n e-foldings of the scale factor. 

As noted above, in an inflationary phase the scale factor accelerates in time, a > 0. The end 
of the inflationary epoch thus occurs at the transition from a > to a < 0. We denote the field 
value at the end of inflation by (p e nd- We find that the value of <f> e nd is virtually insensitive to 
where the field started rolling on the potential, fa. In Paper I and Sect. IIIA, we used fa, the 
value of the field at the end of the SR epoch, as an estimate of the end of inflation. Comparing the 
correct value fa n d with the approximate value fa, we find that the error is only 1% for / ~ Mpi, 
10% for / = O.lMpi, and rapidly gets large for smaller values of /. In particular, no slow rollover 
regime exists for / < Mpi / y/A8ir, and yet for small enough values of fa , significant inflation can 
still occur. In practice, however, the small difference between fa and the exact result fa n d shown 
in Figure 1 is irrelevant, since, as we show below, values of / smaller than 0.3Mp; are excluded 
for other reasons. 

For a given initial value of the field fa (or y\), the solution to Eq. (3.10) tells us the total 
number of inflation e-foldings of the scale factor, N(fa) (where the end of inflation is defined by 
the condition R = 0). Figure 2 shows the number N(fa,f) of e-foldings as a function of the initial 
value of the field fa for different choices of the mass scale /. One can see that, for fa/ f < 1, the 
dependence is almost exactly logarithmic, 

N(fa) = A-B\n(fa/f). (3.11) 
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In the limit of small (f>i/f, the analytic SR estimate of Eq. (3.6) implies this same functional 
dependence and provides values for the constants A and B; in particular, B sr = 167r/ ' 2 /M Pl . 
The numerical values obtained for A and B by solving (3.10) are virtually the same as the SR 
estimates if / is near Mpi and start to differ as / decreases. From Fig. 2, one can read off values 
for y™ ax = (j>i iax / f, the largest initial value of the field that can give rise to N((j)™ ax ) = 60 e- 
foldings of inflation. Again, the numerical results for ^™ aa: are nearly identical to the SR estimates 
(shown in Fig. 1) for values of / near Mpi; they differ by ~ 10% for / = Mpi /10, and deviate 
significantly as / approaches Mpi / \[2Ak from above. 

1) Analytic Solution: Small-angle approximation 

The simple logarithmic behavior of the number of e- foldings N(<p\) indicates that an analytic 
approximation can be found, one which differs from the SR approximation and which is more 
useful for smaller values of f/Mpi. In this region of parameter space, the conditions 

y < 1 and v < 1, (3.12) 

always apply during the inflationary epoch, and the equation of motion (3.10) can be approximated 

by 

where we have made the "small angle" (SA) approximation for the trigonometric functions and 
have neglected higher order terms in y and v. Eqn. (3.13) has solutions of the form 



that is, 



where the constant a is given by 



y(n) = yi e an , (3.14a) 
$ = <f> ie a $ Hdt , (3.146) 



1 r. nil/2 3 3 

* = - 2 [V + H -2 = 2 



1 + Fl 



127T/2J 

Thus, the total number N of e-foldings can be written as 



(3.15) 



N = -ln(cf> 2 /f)--ln( ( f> 1 /f) , (3.16) 
a a 

where 4>2 is the value of the field at the end of inflation. Eqn. (3.16) provides us with an analytic 
solution of the same form as Eqn. (3. 11); note that, here, the constant B of Eqn. (3.11) is given 
by B sa = 1/a, which differs in general from the value B sr predicted by the slow rollover approxi- 
mation. However, for large values of /, such that / S> Mpi /y/12ir, the two approximations agree, 
B sa — > B sr = I6nf 2 /M Pl . Comparison with Fig. 2 shows that, unlike the SR approximation, 
the small angle approximation is also in excellent agreement with the numerical results for small 
values of /. 



C. Constraints: Density Fluctuations, Reheating, Sufficient Inflation 

Having studied the evolution of the homogeneous mode cf)(t) of the scalar field and delineated 
the regions of initial field space for sufficient inflation, we now address other constraints the model 
must satisfy for successful inflation, including density fluctuations and reheating. In particular, 
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these phenomena place tighter constraints on the range of allowed scales / and also limit the 
second mass scale A. Since, in Sec. IIIB, we showed that the SR approximation is accurate for 
the parameter range of interest, we shall rely on it throughout this discussion. 



1) Density Fluctuation Amplitude 

Quantum fluctuations of the inflaton field as it rolls down its potential generate adiabatic 
density perturbations that may lay the groundwork for large-scale structure and leave their imprint 
on the microwave background anisotropy [33-36]. In this context, a convenient measure of the 
perturbation amplitude is given by the gauge-invariant variable (, first studied in ref.[36]. We 
follow ref.fll] in defining the power in (, 

2 V P J HOR 27T <\> 

Here, (bpj p)hor denotes the perturbation amplitude (in uniform Hubble constant gauge) when 
a given wavelength enters the Hubble radius in the radiation- or matter-dominated era, and the 
last expression is to be evaluated when the same comoving wavelength crosses outside the Hubble 
radius during inflation. For scale-invariant perturbations, the amplitude at Hubble-radius-crossing 
is independent of perturbation wavelength. 

To normalize the amplitude of the perturbation spectrum, we assume that the underlying 
density perturbations are traced by the galaxy number density fluctuations up to an overall bias 

factor b g , that is, P^ 2 = P X JJ{ /b g . As inferred from redshift surveys, the variance a 2 gal in galaxy 

counts in spheres of radius 8 h~ l Mpc is about unity (where the Hubble parameter Hq = lOO/i 
km/sec/Mpc). For a scale-invariant spectrum of primordial fluctuations with cold dark matter 
(CDM), this implies [11] 

P c 1/2 ~ (3.18) 

"9 

As we shall see below, we will be interested in cases where the primordial spectrum may deviate 
significantly from scale- invariant, and these cases will be discussed in detail in Sec. IV; here, 
we will use the scale-invariant normalization to get an approximate fix on the scale A. (For 
values of / close to Mpi, this approximation is very accurate.) For the scale-invariant CDM 
model, the recent COBE observation of the microwave background anisotropy [37] roughly implies 

7.7 x 10" 5 < P\ j2 < 1.4 x 10" 4 , or 0.7 < b g < 1.3. 

Using the analytic estimates of Sec. IIIA, the largest amplitude perturbations on observable 
scales are produced 60 e- foldings before the end of inflation, where <p = <fi™ ax , and have amplitude 

p i/ 2 „ A 2 / 9_ /8vr\ 3/2 [1 + cos(^rV/)] 3/2 f3 19) 

< ~ M% 1 2tt V 3 J sin(<^//) ' 1 } 

Applying the COBE constraint above to Eq.(3.19), we find, e.g., 

A = 8.8 x 10 15 - 1.2 x 10 16 GeV for / = M Pl (3.20a) 
A = 1.4 x 10 15 - 2 x 10 15 GeV for / = M Pl /2 . (3.206) 

Thus, to generate the fluctuations responsible for large-scale structure, A should be comparable 
to the GUT scale, and the inflaton mass = A 2 // ~ 10 11 - 10 13 GeV. 

We can use this to determine A as a function of /, shown in Fig.l. For an analytic estimate, 
consider the case / < (3/4) Mpj, for which it is a good approximation to take (f>™ ax / 'it f <C 1. As 
a result, in Eqn.(3.19), we have approximately 

i /2 1.4A 2 / /16^\ 3/2 f f \ 

p < \Y) Urv ' (3 ' 21) 
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Now the last term in this expression is obtained by using Eqn. (3.6) with N(<fi 7 [ iax , fa, f) = 60: 



f 



2 siu I y I ( ' X P 



47T/2 



(3.22) 



Substituting (3.22) on the RHS of (3.21) and using Eqn. (3. 18) we find the value of A(/) in terms 
of the bias parameter: 



M. 



pi 



f 



sin 



92_ 

V 



1/2 



cxp 



87T/2 



(3.23) 



Here, the quantity s'm(fa/2f) is determined by the slow-rollover conditions, Eqns. (3.4-3.5) and 
is generally of order unity. The dominant factor in (3.23) is the exponential dependence on f 2 , 
which is responsible for the rapid downturn as / begins to drop significantly below Mpi in the 
curve for A(/) in Fig.l. For completeness, we note that the value in Eqn. (3.23) is strictly only an 
upper bound on the scale A, since the perturbations responsible for large-scale structure could be 
formed by some other (non-inflationary) mechanism. 



2) Density Fluctuation Spectrum 

Using the approximations above, we can investigate the wavelength dependence of the per- 
turbation amplitude at Hubble-radius-crossing and in particular study how it deviates from the 
scale-invariant spectrum usually associated with inflation. Here we give a quick derivation of the 
spectrum, and defer a fuller discussion to Sec. IV. 

Let k denote the comoving wavenumber of a fluctuation. The comoving lengthscale of the 
fluctuation, k , crosses outside the comoving Hubble radius [Ha] -1 during inflation at the time 
when the rolling scalar field has the value fa. This occurs Ni(k) = N(fa, fa, /) e-folds before 
the end of inflation, where N(fa,fa,f) is given by Eqn. (3.6) with fa replaced by fa. The 
corresponding comoving lengthscale (expressed in current units) is 

k~ l ~ (3000/ l - 1 Mpc)exp(A^ / (fc) - 60) . (3.24) 

For scales of physical interest for large-scale structure, Nj(k) > 50; for / < (3/4) Mpi, these scales 
satisfy fa/ f <C 1. In this limit, comparing two different field values fa x and 4>k 2 , from Eqn. (3.6) 
we have 

fa 2 c fel exp (- , (3.25) 

where ANj = N I (k 2 ) — Nj(ki). Thus, using Eqns. (3. 19) and (3.21), we can compare the pertur- 
bation amplitude at the two field values, 



1/2^ 

~ exp ( - ~ x,1 *"Pi ) . (3.26) 



< 2 ) fc2 fci V 16*/ 



Now, from Eqn. (3. 24), we have the relation ANj = ln(/ci//c2) (more precisely, ANj = ln(/ci ^2/^2^1) 
and we approximate Hi ~ H 2 ); substituting this relation into (3.26), we find how the perturbation 
amplitude at Hubble radius crossing scales with comoving wavelength, 



( S -P) ~ (Pl /2 ) k ~ k -^/^f 2 (3.27) 

V P J HOR.k 
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By comparison, for a scale-invariant spectrum, the Hubble radius amplitude would be independent 
of the perturbation lengthscale A; -1 ; the positive exponent in Eqn.(3.27) indicates that the PNGB 
models have more relative power on large scales than scale-invariant fluctuations. 

It is useful to transcribe this result in terms of the power spectrum of the primordial per- 
turbations at fixed time (rather than at Hubble-radius crossing). Defining the Fourier transform 
8k of the density field, from Eqn.(3.27) the power spectrum is a power law in the wavenumber k, 
(\8k\ 2 ) ~ k n % where the index n s is given by 

M 2 

For comparison, the scale-invariant Harrison-ZePdovich-Peebles-Yu spectrum corresponds to n s = 
1. For values of / close to Mpi, the spectrum is close to scale- invariant, as expected; however, as 
/ decreases, the spectrum deviates significantly from scale-invariance-e.g., for / = Mpij \[&k = 
0.2Mpi, the perturbations have a white noise spectrum, n s = 0. In sec. IV, we explore the 
implications of models with power law primordial spectra in depth. 



3) Quantum Fluctuations 

For the semi-classical treatment of the scalar field used so far to be valid, the initial value 
of the field should be larger than the characteristic amplitude of quantum fluctuations in (ft, i.e., 
(fti > A(ft = H/2-k. In particular, requiring that quantum fluctuations do not reduce the number 
of inflation e- folds below 60 implies that the condition (ft™ ax > H/2ir must be satisfied. Using the 
SR. approximation and Eqn.(3.22), we find 

H/2n ^ J_ 1/2 /MpA 2 _ 10^ /MpA 2 

W°* ~ 3?r C \ f J ~ b g \ f ) ' (3 ' 29) 



Since this ratio is very small over the parameter range of interest, this constraint places no 
significant restrictions on the model. For example, this constraint requires that (fti/f > 10~ 7 for 
/ = Mpi and (fti/f > 6 x 10~ 9 for / = Mpi/2, while the corresponding values of (ft 7 [ tax / f are 0.63 
and 9.4 x 10~ 3 . Even if (ft\ is at some stage smaller than this constraint, one would expect that 
quantum fluctuations might eventually bring the field into the semiclassical regime, so inflation 
would begin, if it was sufficiently spatially coherent. 



4) Probability of Sufficient Inflation 

Armed with the numerical and analytic results above, we now calculate the a posteriori 
probability of sufficient inflation. We consider the universe at the end of inflation, and calculate 
the fraction V of the volume of the universe at that time which had inflated by at least 60 
e-foldings: 

JJrL* #iexp[3iV(0i)] 

V = l- . (3.30) 

Jl 1 ^ # ie xp[3iV(<M] 

Here, the lower limit of integration in the denominator is the limit of validity of the semiclassical 
treatment of the scalar field; the initial value of (ft must exceed its quantum fluctuations, (ft\ > 
A(ft = H/2tt. We will use the form for N((ft 1 ) given by Eq. (3.11) to evaluate the integrals 
appearing in Eq. (3.30). As shown previously, this approximate form for N((ft\) is only valid 
for (fti/f < 1. However, we will assume that it holds over the entire range of integration; in the 
Appendix, we show that the resulting errors are small. Our basic result is that the a posteriori 
probability for inflation is essentially unity for / larger than the critical value f c ~ 0.06Mp;. As 
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/ drops below this value, the probability given by Eq. (3.30) rapidly approaches 0. To illustrate 
this result, we evaluate the integrals in (3.30); both are of the form 



dyie 



3A -3B\n Vl 



;l , (3-31) 

~ 35^11 V ~ [ n> /' 



where e is a small number, e = H/2irf or y™ ax , and B is the /-dependent coefficient appearing 
in Eqn.(3.11). If 3B > 1, the integral / is dominated by the lower end of the range of integration 
and only the first term in Eq. (3.31) is significant. In this case, the probability V is given by 

/ H \3B-1 

V = l-( ) , (3.32) 

where 3-B — 1 is positive. The probability of sufficient inflation is close to unity as long as the 
ratio in brackets H/2ir(ft r [ iax is small; however, this is guaranteed by Eqn.(3.29). Combining Eqs. 
(3.29) and (3.32) yields the probability: 



V > 1 



10" 5 ( M 



L b g 



i pi 



2 



3B-1 



(3.33) 



This expression is valid provided that 35 — 1 is positive and not extremely close to 0. 

As the value of 3B decreases toward unity, the probability V decreases and the approximation 
leading to (3.32) begins to break down. As a reference point, consider the special case 3B = 1; then 
the integral / = e 3j4 ln(-7r/e) and the probability V ~ 0.05. As B decreases further, the integral 
in Eq. (3.31) obtains most of its contribution from the upper end of the range of integration 
and hence both integrals appearing in Eq. (3.30) have nearly the same value. As a result, the 
probability V rapidly approaches 0. 

To summarize, we find that the probability V of sufficient inflation depends primarily on 
the value of the coefficient B appearing in Eq. (3.11), which in turn determines the number of 
e-foldings of the universe as a function of the initial value <p\ of the field. For B > 1/3, the 
probability V is nearly unity; for B < 1/3, the probability V quickly approaches 0. In the SR 
approximation, B k, IQnf ' 2 /Mp t , which would imply a critical value j & c r = j(B = 1/3) = l/\/487r. 
On the other hand, the numerical calculations yield the critical value of the mass scale f c = 0.058. 
This discrepancy is traced to the fact that the SR approximation is invalid for such small values 
of /. In this case, the "small angle" approximation discussed in Sec. III.B is more appropriate; 
using Eqn.(3.15), we can analytically determine the critical value of / for which B sa = 1/a = 1/3, 

This is in excellent agreement with the value found numerically. 



5) Reheating 

At the end of the slow-rolling regime, the field (f> oscillates about the minimum of the potential, 
and gives rise to particle and entropy production. The decay of (j) into fermions and gauge bosons 
reheats the universe to a temperature Tjjh — (45/4^r 3 g*) 1 / 4 ^f^Mp~l, where g* is the number of 
relativistic degrees of freedom. On dimensional grounds, the decay rate is F ~ g 2 m ( f > 3 /f 2 = 
g 2 A 6 /f 5 , where g is an effective coupling constant. (For example, in the axion model [6,7], 
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g oc aEM for two-photon decay, and g 2 oc (m^/m^) 2 for decays to light fermions ip.) Thus, the 
reheat temperature is 

/ 45 \ 1/4 5 A 3 (M Pl \ l/2 

For example, for / = M P i, using (3.20a) for A, and taking g* = 10 3 , we find T RH ~ 10 8 g GeV, 
too low for conventional GUT baryogenesis, but high enough if baryogenesis takes place through 
sphaleron-mediated processes at the electroweak scale. Alternatively, the baryon asymmetry can 
be produced directly during reheating through baryon- violating decays of (ft or its decay products. 
The resulting baryon-to-entropy ratio is ris/s ~ eT^n/m^ ~ egA/f ~ W~ 4 eg, where e is the 
CP-violating parameter; provided eg > 10~ 6 , the observed asymmetry can be generated in this 
way. 

We saw above that the amplitude of density perturbations produced during inflation yields 
a bound on the scale A as a function of the fundamental scale /, eqn.(3.23). We can use this to 
express Trh as a function of / (which depends only weakly on g and g*); requiring that this be 
sufficiently high for some form of baryogenesis leads to an important lower bound on the scale /, 
which as we shall see below, is more restrictive than the a posteriori bound above and comparably 
restrictive with the microwave anisotropy bound on the perturbation spectrum to be discussed in 
Sec. IV. Since we will be interested in a lower bound on /, we consider the case / < (3/4)Mp; so 
that eqn.(3.23) applies. Substituting (3.23) into (3.35), we find the reheat temperature 



10'° GeV flOB\ l ' , fMp l \ i . 3/2 ( <h , 
T R „ = — ^ (-) [—) s,n/ (^-|exp 



'g 



45M| 
Svr/2 



(3.36) 



The important point here is that the reheat temperature drops exponentially as / drops well below 
Mpi. For baryogenesis to take place after inflation, at a minimum we should require Trh > 100 
GeV, the electroweak scale. From eqn.(3.36), this leads to the lower bound 

f 

> 0.28 . (3.37) 
Mpi ~ 



(Here, we have set g = 1 and g* = 100, but this limit depends only logarithmically on g and 
g*.) In terms of the density perturbation spectrum given in Eqn.(3.28), if inflation produces the 
dominant fluctuations on all scales, then this reheating constraint implies n s > 0.5. 

One additional point concerning reheating in these models deserves mention. In the string 
models of §11, the axion couples predominantly to the hidden sector; in such inflation models, one 
might then worry that reheating would take place more efficiently in the hidden as opposed to 
the ordinary sector. (This would not be a concern in models without a hidden sector, such as 
those patterned after technicolor.) In practice, this is not an insurmountable obstacle for these 
models, because gravitational interactions lead to an effective coupling between the hidden sector 
inflaton and the ordinary sector particles. Furthermore, for / ~ Mpi, the gravitationally induced 
decay rate to ordinary particles, T ~ m^/M Pl , is comparable to the axion's decay rate to the 
hidden sector. Thus, we would expect the two sectors to reheat to comparable temperatures. It 
is then easy to imagine a subsequent entropy-producing ordinary particle decay which heats the 
ordinary sector relative to the hidden sector, so that the contribution of the hidden sector to the 
total energy density at the time of big bang nucleosynthesis is negligible. 



6) Initial Spatial Gradients 

In the previous discussion, we have focussed on the evolution of a nearly homogeneous scalar 
field 4>(t). However, since we expect the field initially to be laid down at random on scales larger 
than the Hubble radius, spatial 'Kibble' [30] gradients will be present on these scales. For inflation 
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to occur, it is necessary that the stress energy tensor averaged over a Hubble volume be dominated 
by the potential V(4>), not by gradient terms ((di<p) 2 ). (This is of course a concern for all models 
of inflation, not just those considered here.) In paper I, we addressed this issue at some length, 
and argued that, when the universe has cooled to the temperature T ~ A at which inflation would 
otherwise begin, the energy density contributed by field gradients would be at most comparable to 
that in the potential. (During the prior radiation-dominated epoch, the gradient energy density 
scales like radiation, p gra d ~ (di<fr) 2 ~ f 2 /t 2 ~ T 4 , where the last equality assumes / ~ Mpj; 
thus, at T ~ A, we expect p gra d ~ A 4 ~ V(4>).) Since these gradients rapidly redshift away with 
the subsequent expansion, they would typically delay only slightly the onset of inflation. 

Here, we point out that the canonical PNGB model has an additional automatic feature which 
can ensure that spatial gradients in the PNGB field are negligible at the onset of natural inflation. 
Namely, if (f> is the angular component of a complex field <]?, as in the model of Eqn.(2.3), then 
the heavier, radial component of $ can generate an earlier period of inflation as it rolls down its 
potential. If the later angular inflation leads to more than 60 e-folds of growth in the scale factor 
(as we have been assuming), then the only important effect of the earlier inflation epoch would be 
to rapidly stretch out spatial gradients in the angular <j) field. (This point was stressed to us by A. 
Linde, private communication.) Furthermore, as we show below, the earlier inflation period does 
not require another small coupling constant. In particular, for the model of Eqn.(2.3), for a broad 
range of initial conditions, radial inflation takes place even if the complex scalar self-coupling A is 
of order unity. In addition, only a small number of radial inflation e-folds is required to efficiently 
damp spatial gradients in <p. 

In the usual way, we can decompose the complex field $ into two real radial and angular 
components rj and <p, 

$ = e ^//J? (3.38) 
v2 

Consider the evolution of the radial mode 77 in the potential (2.3), V(rj) = (A/4)(r/ 2 — f 2 ) 2 (in 
general the radial and angular motions are coupled; however, since the radial mode is much 
heavier, its evolution can be approximately decoupled). Analyzing this motion in a manner 
analagous to §111. A, and using the fact that / is comparable to Mpi, we see that some amount 
of radial inflation is expected provided the initial value of r\ is sufficiently far from its minimum 
(r/) = f. In fact, this initial period of inflation will be generic as long as gradient terms in the 
r\ energy density do not dominate over the potential V(rj) near the Planck scale and the initial 
value of i] is not very close to /. For example, for / = Mpi, if the initial value r/i of the radial 
field is greater than 2Mpi, then in rolling to its minimum it will generate at least 5 e- foldings of 
'chaotic' inflation, and angular gradients would be stretched by a large factor. Alternatively, if 
r/i < 0.3Mpi, the universe would experience about the same number of e- foldings of 'new' inflation 
as the field rolls from near the local maximum of the Mexican hat at the origin. We note that, 
for a potential of the form (2.3), for / near Mpi the SR condition holds over some range of rj, 
independent of the value of the coupling A (just as Eqn.(3.6) does not depend on A). Therefore, 
radial inflation takes place even if A is large. The density fluctuations produced during this phase 
are on unobservably large scales if the subsequent angular inflation lasts for at least 65 e-folds 
of expansion, so there are no strong constraints on A arising from density fluctuations and the 
microwave anisotropy. One should, however, require \/A/£ < 1 to avoid fluctuations of order 
unity on the Hubble radius, since these would pinch off into black holes. 



IV. POWER LAW SPECTRA AND LARGE-SCALE STRUCTURE 



Recent observations of large scale galaxy clustering and flows suggest that there is more 
power on large scales than the n s = 1 scale invariant spectrum gives for 'standard' cold dark 
matter dominated universes (CDM models). In this section, we show the degree to which varying 
the index n s , where the primordial power spectrum \Sk\ 2 ~ k ns , while keeping all other features of 
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the CDM model fixed, helps solve this large scale structure dilemma. We have shown that natural 
inflation will generate such a power law perturbation spectrum over a wide range of wavenumbers, 
in particular over the waveband that we directly probe with observations of large scale galaxy 
clustering and of microwave background anisotropies. We demonstrate this in more detail in Sec. 
IV. A below. In addtion, other inflation models, such as 8those with exponential potentials and 
many versions of extended inflation, also predict power law spectra which can deviate from scale 
invariant. In Sec. IV. B, we show that current data on microwave anisotropies and large-scale flows, 
and the requirement that structure forms sufficiently early, constrain n s to be > 0.6 for CDM 
models, but values < 0.6 are needed to explain the large scale clustering of galaxies. The reason 
we put the CDM model under such scrutiny rather than other inflation-inspired models, apart 
from its having dominated the theoretical scene for the past decade, is that it is a minimal model, 
in the sense that it requires the least number of assumptions to specify it. For the 'standard' CDM 
model, one assumes a flat geometry for the Universe with Q ~ 1 in non-relativistic particles and 
takes h~ 0.5, where h is the Hubble constant Hq in units of 100 km s _1 Mpc _1 . (For values of h 
larger than this, if f2 = 1 the Universe would be younger than the inferred ages of globular cluster 
stars.) We assume negligible baryon abundance, Qb *C ^, in the following; a small value of 
is indicated by primordial nucleosynthesis constraints (< 0.07). The rest of the non-relativistic 
matter is in cold dark matter relics, Q c dm = ^ — Qb- Since the large scale structure dilemma has 
been with us in one guise or another since the early 1980s, a major line of research over the past 
decade has been to invent models with scale invariant primordial spectra that have more power 
than the n s = 1 CDM model does on large scales. These 'nonstandard' n s = 1 models include 
models with a non-zero cosmological constant, a larger baryon density £Ib than that inferred from 
standard nucleosynthesis, and mixtures of hot and cold dark matter, to name just a few. Often 
somewhat baroque from the particle physics prespective, such alterations would all result in more 
stringent constraints on n s if we allow it to vary than the ones we derive for the CDM model. 
(Indeed there are models that require the effective n s to be 3> 1, such as the isocurvature baryon 
model, but this is certainly not an outcome of natural inflation.) 

A. Inflation Models and Power Law Spectra 

Before turning to the data, we first show explicitly how tiny the deviations from a power 
law form are for natural inflation, and that Eqn.(3.28) for n s is highly accurate. We also discuss 
the form that n s takes for other popular models of inflation such as power law, extended, and 
chaotic inflation. Since we are dealing with spectra that can change somewhat with wavenumber, 
we define a 'local' (i.e., fc-dependent) spectral index n s [k) by 



where the C-power spectrum P^(k) introduced in Sec. Ill provides a better measure of the post- 
inflation spectrum than does the density power spectrum. The quantity £, the variation of the 
3-space volume on uniform Hubble parameter hypersurfaces, is gauge- and hypersurface-invariant, 
whereas the density is neither. 

1) Natural Inflation 

For natural inflation this local index is 



Here Nj(k) is the number of e- foldings between the time when the inverse wavenumber k 1 first 
exceeded the comoving Hubble length (the first 'horizon crossing') and the end of inflation. For 



n s (k) = 1 + dlnP c (fc)/dlnlfe , 



(4.1) 
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(4.2) 
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waves on scales of observable interest, Ni(k) ~ 50 — 60, so the fator in large brackets is always 
very close to unity over the entire range of values of / we are considering. 

The derivation of (4.2) is very similar to that given in Sec. Ill, so we just sketch the steps 
here. From Eqns.(3.17) and (4.1), we must evaluate n s — 1 = 2dln((3 /2ir) H 2 /\cp\)/ din Ha (since 
k = Ha at horizon crossing). If we use the slow roll approximation for and H, we have 



n s (k) 
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(4.3) 



Here 4>k is the value of the scalar field at which k = Ha. As in Sec. Ill, we have taken the positive 
sign for the potential (1.1). The 'end of inflation' occurs when the scalar field kinetic energy grows 
to the value 4> 2 = V, i.e., the Universe passes from acceleration to deceleration (cf. Eqn.3.3). At 
this point, the expansion rate H ~ (3/2) 1 / 2 H sr , where the Hubble parameter during the SR 
epoch, H sr = y/8irV(4>)/3Mp l . As a result, the end of inflation can also be expressed as the time 
when \4>\ = HMpij\f^K. Approximating cf> by the slow roll result, 4> = (— Mp 2 /4ir)dH sr /d4>, 
one finds that inflation ends when the field reaches the value <ft en d — 2/ arctan [\/247r// Mpi ] . 

■ 2 

(In Sec. Ill, we defined the end of SR to occur when <fi 2 = (f>2 = 2V, which gave the factor 
\/487r in the argument of the arctan [Eqn.(3.5)] rather then the \/2Att found here. The numerical 
computations of <\> en d discussed in Sec. III.B are best fit by a \/347r factor in the argument, so 
eq.(3.5) or the value given here give about the same accuracy. In any case, in Eqn.(4.2) this factor 
is multiplied by the exponential suppression factor exp [-(Mp/ 2 /87r/ 2 ) Nj(k)].) By solving the 
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This expression generalizes Eqn.(3.7); when it is substituted into Eq.(4.3), Eq.(4.2) is obtained. 

Defining k en d to be the wavenumber that equals (Ha) en d at the end of inflation, and using 
the fact that Nj(k) = ln(H (cf)k)k en d / H (4> en d)k) , the relation between Ni(k) and k is given by 
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Thus between the current Hubble length /c _1 ~ 3000h _1 Mpc and the galactic structure length 
scale, k~ l ~ 0.5 h _1 Mpc, the range which encompasses all of the large scale structure observations, 
Nj(k) only changes by about 10. Since Nj(k) only enters the exponentially suppressed terms in 
(4.2), the index n s is quite constant at 1 — Mp 2 /(87r/ 2 ) over observable scales. 



2) Exponential Potential Inflation 

Although we view natural inflation as the best motivated model for obtaining power law 
indices below unity, other possibilities for getting n s {k) significantly different from unity have 
been widely discussed in the literature. Power law inflation 8[38,39] (in which the scale factor 
grows as a large power p of the time, a oc t p , instead of quasi-exponentially) is the simplest 
example of a model which predicts power law spectra. It is realized with an exponential potential 
of form V = Voexp[— ^\Qir/p cf)/Mpi], and has 

n s = 1 - . (4.4) 

p — 1 
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The deceleration parameter of the universe, q = —ad/a 2 is q = — (1 — for power law inflation. 
In order to have a viable model of inflation, the universe must pass from acceleration, q < 0, to 
deceleration, q > 0, so that it can reheat, hence it is essential that p evolves, with inflation ending 
when p falls below unity. Thus, although power law inflation models are instructive since they are 
analytically simple, the exponential part of the potential can only be valid over a limited range of 
the evolution. Indeed, it is often convenient to characterize potentials that are not exponentials 
by an index p defined by y/Airp = HMpi/\4>\, which reduces to the p in the exponential potential 
for that case. However, in these models structure on observable scales may be generated in a 
regime where p varies with k rather than being constant. Even so, power law approximations are 
often locally valid, even when rather drastic potential surfaces are adopted to 'design' spectra. 
Some examples of cases where n s changes considerably over the observable window of large scale 
structure are given, for example, in [40,11]. 



3) Extended Inflation 

Extended inflation also leads to a power law form over a wide band in fc-space [41]. In ex- 
tended inflation, a Brans-Dicke field, whose inverse is an effective Newton gravitational 'constant', 
is introduced as well as an inflaton. The analysis of [41] showed that the power law index can be 
simply expressed in terms of the Brans-Dicke parameter to (the coefficient of the kinetic term of 
the Brans-Dicke field), 

As far as density fluctuations are concerned, the model just mimics a power law inflation one. 
Indeed, the fluctuation spectrum is most easily computed in a conformally-transformed reference 
frame in which the log of the Brans-Dicke field experiences an exponential potential with p as 
given in Eq.(4.5) [41], yielding the n s relation through Eq.(4.4). We note that, in most versions of 
the theory, a value of to < 18 — 25 is needed to avoid an excessive CMBR anisotropy due to large 
bubbles, which implies that the spectrum deviates from scale invariant, n s < 0.77 — 0.84. At the 
same time, it is also necessary that the effective value of uj must have evolved to a high number 
(> 500) by now in order to satisfy solar system tests. This can be arranged by, e.g., giving the 
Brans-Dicke field a mass or by other means, but at the cost of complicating the model. 



4) Chaotic Inflation 

Although references [40,11] probed how dramatic the breaking of scale invariance could be in 
terms of the fluctuation spectra over our observable waveband, the main conclusion was that plau- 
sible models of inflation were much more likely to lead to quite smooth breaking over the observ- 
able range. We illustrate the level of breaking of scale invariance expected for the popular chaotic 
inflation models. We assume power law potentials of the form V(4>) = X e Mpi 4 ((f>/Mpi) 21 ' /(2is), 
where the power v is usually taken to be 1 or 2. A characteristic of such potentials is that the 
range of values of <fi which correspond to all of the large scale structure that we observe is actually 
remarkably small. For example, for v = 2, the region of the potential curve that corresponds to 
all of the structure between the scale of galaxies and the scales up to our current Hubble length 
is just 4Mp; < 4> ~ 4AMpi [11]. Consequently, the Hubble parameter does not evolve by a 
large factor over the large scale structure region and we therefore expect near scale invariance. 
Although this is usually quoted in the form of a logarithmic correction to the (^-spectrum, a power 
law approximation is quite accurate. Following exactly the same prescription used to evaluate 
Eq.(4.2), we have 

For waves the size of our current Hubble length we have the 8familiar Ni(k) ~ 60, hence n s ps 0.95 
for v = 2 and n s ps 0.97 for v = 1 (massive scalar field case). The relation between Nj(k) and k 
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is given by 



In 




(4.7) 



where k en( i is the wavenumber that equals Ha at the end of inflation. Thus, over the range from 
our Hubble radius down to the galaxy scale, n s decreases by only about 0.01. 

B. Implications for Large-Scale Structure 

We have discussed various inflationary models (natural, power law, extended, and chaotic), 
which give rise to density perturbation spectra of the form \5k\ 2 ~ k ns , where n s < 1. We now 
turn to their implications for large scale structure. 

1) Galaxy and Cluster Clustering 

Although the amplitude of the fluctuations is known once all aspects of the inflaton potential 
are specified, it is more convenient to normalize the spectrum to the level of clustering we observe 
and use that to restrict particle physics parameter ranges, as in Sec. III. We normalize the ampli- 
tude of the density perturbation spectra by setting the rms fluctuation in the mass distribution 

within spheres of radius 8h _1 Mpc, a p ,g, = ((5M / M) 2 ) 1 ^ 2 _ sh _ 1Mpc , to be <7g. The rms fluctuation 

in galaxy counts on this scale in the CfA survey is unity. The quantity fog = a^ 1 is sometimes 
called the 'biasing' factor, since roughly if fog ~ 1 we expect that galaxies would be clustered like 
the mass distribution while if fog > 1 galaxies would be more strongly clustered than the mass; 
this point is discussed in more detail below. For standard CDM models with n s = 1, <7g was 
thought to lie in the range 0.4 — 1.2 before the recent COBE measurement. 

The evolved power spectra of the linear CDM density fluctuations, da 2 p /d In k = k 3 (\5k{to)\ 2 ) /2tt 2 , 
are shown in Figure 3(a), for spectral indices ranging from n s = —1 to 1. A transfer function T(k) 
relates the primordial spectrum \5k(ti)\ 2 oc k ns to the present spectrum, |5fc(to)| 2 = T (k)\Sk(ti)\ . 
For the CDM transfer function, we use the fitting formula given in Appendix G of BBKS [42], 
which is highly accurate in the — > limit, but is somewhat modified for the ~ 0.05 values 
more appropriate from nucleosynthesis. The spectra are in units of <7g. The spectra are plotted 
in this way to provide a measure of the contribution of a band around the given wavenumber to 
the overall rms density fluctuations; the ordinate roughly gives (1 + z n i(k))/as, where z n i(k) is 
the redshift at which the rms fluctuations in the band become nonlinear. Notice that there is a 
peak in the CDM spectrum for n s < 1. This indicates that there is a characteristic scale, roughly 
the peak, associated with the first objects that form [43]. A potential problem with these models 
that is immediately apparent from Fig. 3(a) is that the redshift of galaxy formation is lower than 
that for the scale invariant model, which, for small n s , can lead to grave difficulties in explaining 
why there are quasars at z ~ 5. We discuss this point more fully below. 

To relate such a linear density perturbation spectrum to galaxy clustering, one must generally 
do iV-body calculations. However, on large scales, the waves evolve in an essentially linear fashion, 
and there is an excellent approximation which relates the power spectra of galaxies and clusters 
of galaxies (if they arise from any function of the Gaussian process through which perturbations 
arose) to that of the density field. This relation is an extension [44] of the theory which identifies 
galaxies and clusters with appropriately selected peaks of the initial density field [42, 45]. For 
scales large compared with the local processes that define these objects and large enough that the 
waves are evolving in the linear regime, the power spectra for galaxies and clusters are linearly 
proportional to the density power spectrum, with the proportionality constants defining 'biasing 
factors', b 9 for galaxies and fo c for clusters (Cf. Sec. III.C.l): 



da 2 



b 2 p - c 

g dlnk ' din A; 




(4.8) 



dink 
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In Figure 3(b), which focuses on the region of fc-space in 3(a) probed by large scale structure 
observations, we compare the predicted galaxy spectrum with large scale clustering data from the 
QDOT and UC Berkeley IRAS surveys and (less directly) the APM survey. The spectrum is in 
units of bgCTg- In the conventional BBKS peaks approach to biasing [42], we would have b g = 1/<t 8 , 
which is why a^ 1 , the inverse of an amplitude measure, is often referred to as a biasing factor 
(e.g., in Sec. Ill) . In general, b„ will differ from galaxy type to galaxy type and there is no 
clear reason why we should suppose that b g = a^ 1 ; nonetheless, it is rather remarkable that this 
assumption appears to give the correct amplitude for galaxy clustering. However, we note that 
the slight differences in the power spectrum levels for the 3 surveys could be simply explained 
with slightly differing b g 's. To compare with the data in the nonlinear regime of the spectrum, 
< 5cr 8 h _1 Mpc, iV-body computations are needed. However, just from the linear regime it 
would appear that spectral indices in the range 0-0.6 are much preferred over the scale invariant 
value of unity. (This point would appear even more dramatic had we forced the models to agree 
with the data at the 8 h~ x Mpc normalization scale.) 

Probably the most reliable indication of excess large scale power is the angular correlation 
function of galaxies, w gg (9), inferred from deep photometric surveys. Although the angular corre- 
lation function suffers from having only two- rather than three-dimensional information, it gains 
enormously since angular surveys currently involve a few million galaxies, while three-dimensional 
(redshift) surveys are still limited to samples of several thousand galaxies. Two groups have now 
independently catalogued the galaxies of the Southern Sky and have derived w gg (9ys in agree- 
ment with each other. A Northern Sky survey is also in basic agreement. Bond and Couchman 
[44] showed that the angular correlation function at large angles can be evaluated using the lin- 
ear power spectrum for galaxies, although nonlinear effects substantially modify the estimates at 
small angles; they also showed how to evaluate the correlation function directly from the power 
spectrum. We applied these techniques to the power spectra of Figure 3(b) to compare w gg as we 
vary n s with the APM results in Figure 4. The dots denote the APM data for various magnitude 
intervals, scaled back to the depth of the Lick survey [46]. The spread is considered to provide a 
rough indicator of the error level. Although there is a certain amount of vertical freedom in fitting 
the theory to the data, from the overall scale b g as, it is clear that < n s < 0.4 is required if we 
are to take the spread of dots as an error estimate. It was this graph that led to the conclusion 
given in Bond [47] that this was the allowed range. However, estimates for various corrections to 
the APM catalogue such as those from plate errors and variable absorption by Galactic dust may 
revise w gg {6) downward slightly, and the hatched region is now expected to be allowed by the 
data [46]. Thus, for this paper, we consider the allowed range to be < n s < 0.6. We note that 
this fit has been done with a CDM spectrum with h=0.5 and VLb ~ 0. If we can contemplate h 
as low as 0.4 or Q,b as large as 0.1, then n s « 0.7 is feasible as well. 

The high degree of clustering of clusters has been a puzzle since the early 1980's. The 
correlation function of rich clusters was thought to be enhanced by a factor of about 11-16 over 
the level of galaxy clustering, assuming both have the same power law behaviours [48]. The sample 
from which most of the estimates of clustering were derived was the Abell catalogue, which has 
been criticized on a number of grounds. The main problem seems to be the projection effect, in 
which clusters at different redshifts superimpose upon one another, leading one to believe that the 
clusters are more massive than they truly are. Recently two redshift surveys of clusters identified 
using the Southern Sky galaxy surveys estimate correlations about half as large as the original 
values, and have shown that they are not subject to contamination by projection effects. These 
new values are roughly compatible with the levels expected if one uses the power spectra suggested 
by the galaxy clustering data [49]. Provided we are in the linear regime, Eq.(4.8) shows that the 
correlations should be in the ratio (b c /b g ) 2 . A rough estimate for this ratio can be obtained using 
the methods of [45] for a peak model of clusters, in which one can determine the combination 
(b c — l)as just from the abundance of clusters; it is about 2.1. Thus, (b c /b g ) 2 ~ (2.1+ag) 2 / (b g as) 2 ■ 
Taking b g = a^ 1 , the enhancement factor ranges from 6 to 10 as <jg ranges from 0.5 to 1. Thus, if 
the new cluster correlation functions prove to be valid, they can also be explained with the same 
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range of n s as the w gg data indicates. 



2) Constraints from Microwave Background Anisotropies 

We now determine the range of as as a function of n s suggested by the COBE observations of 
microwave background anisotropy with the Differential Microwave Radiometer experiment [37]. 
The DMR team have given data for the rms fluctuations on the scale of 10°, <rr(10°), the sum of 
the squares of the components of the quadrupole moment tensor, which we denote here by a\ e=2 , 
and estimates of the correlation function with the dipole and quadrupole contributions removed. 
Here, we express all of these values in units of AT/T, by dividing their results by the background 
temperature, 2.736K. 

The fwhm of the DMR beam (7°) is sufficiently large that it is quite accurate to assume for 
the adiabatic fluctuations of interest here that the microwave background anisotropies arise from 
curvature fluctuations experienced by the photons as they propagate through photon decoupling 
to the present (Sachs- Wolfe effect). If we assume that the universe is matter dominated from 
photon decoupling to the present, the variance Cg of the multipole coefficient ag m in the spherical 
harmonic expansion of the radiation pattern (see e.g., ref. [50]), is given by 

Cg = (M 2 ) = ^ dink ^ i|(fcro) , (4.9a) 

where jg is a spherical Bessel function and tq is the comoving distance to the photon decoupling 
region, To — 2i7 ~ 1 ~ 6000h _1 Mpc. The comoving wavenumber k is referred to current length 
units. The gravitational potential spectrum is related to that for the density by 

da%/dlnk = ({3/2)H 2 k- 2 ) 2 da 2 p /dlnk . (4.96) 

Although we used Eq.(4.9a) directly to evaluate the temperature power spectrum Cg, for power 
law spectra on the large scales that COBE probes, there is a simple expression in terms of Gamma 
functions [50] and the quadrupole power C2: 

T[i+ (" 3 - 1 ) ipr (9-"8) i 
Cl = C \[i + ^A ] ^A ] ' (410) 

for i>2. In terms of Cg, the rms value expected in each multipole for COBE is 

<4i = ^Cg T 2 , (4.11) 

where Tg is a filter appropriate to their beam, and is approximated by a Gaussian 
Tg = exp[-0.5(£ + 0.5) 2 /(£ dmr + 0.5) 2 ] , l dmr w 19 , 

where ldmr corresponds to 7° fwhm. 

The strongest result to use for estimating the amplitude as is provided by <tt(10°), which 
the COBE team determined by evaluating the intrinsic sky dispersion after further smoothing 
their data with a 7° fwhm Gaussian filter. To compare with this, we calculate the average value 
that our theoretical model predicts for this, 



4(10°) = £ T 2 ga 2 Tl . (4.12) 
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The extra filtering by Tf brings the total smoothing up to a total of 10°. Since the realization of 
the Universe that we observe involves a specific set of multipole coefficients drawn from (Gaussian) 
distributions with variance Cg, there will be a theoretical dispersion in the values of crf^lO ), what 
the COBE team refers to as cosmic variance. For 0^(10°), we have 

<[A4(10°)] 2 > =2^^ T [^ 2 4,] 2 . (4.13) 

An excellent fit to our calculation of Eqs.(4.12,4.13) is 

a T (10°) = 0.93 x 10" 5 a 8 e 2 ' 63(1 -" s) [1 ± 0.1e a42(1 - ns) ] . (4.14) 

(Since the error is for the square, <7j.(10°), there is a slight asymmetry between the upper and 
lower error bars for 0^(10°) which we have included in Figure 5.) Eq.(4.14) is to be compared 
with the DMR result, including their '1 sigma' errors, 

{(TT(W°)] dmr = 1.085 x 10 -5 [1 ± 0.169] . (4.15) 

(These errors should be slightly enhanced since the detected large scale anisotropy can lead to 
bigger fluctuations in <tt(10°) than one would get solely using single pixel errors, as the DMR 
team did. This appears to be a sufficiently small correction that it can be ignored.) The combined 
error is therefore about 20% for n s = 1, rising slightly for lower values, hence 



^8 



= 1.17 e -2.63(l-n s ) [ 1±0 .2] . (4.16) 



The allowed region for ag as a function of n s using our computed values, is shown in Figure 5. In 
particular, for n s < 0.6, the DMR result requires a 8 < 0.5. 

However, we caution that this value is for the Qb = limit. With the value Qb ~ 0.06 
favoured by primordial nucleosynthesis, <jt(10°) rises by about 15% and ag drops by this amount. 

The quadrupole determination by the DMR team is not nearly as restrictive, because the 
'cosmic variances' as well as the DMR error bars are quite large. Integrating Eq.(4.9a) over all 
k > 10" 4 h _1 Mpc for C 2 , we obtain 

a T ,e=2 = 0.46 x 10" 5 a 8 e 2 - 94 ^""^ [1 ± 0.3] , (4.17) 

to be compared with 

[(T T ,e=2}dmr = 0.475 x 10" 5 [1 ± 0.31] , (4.18) 

hence 

ag « 1.02e- 2 - 94(1 -" s) [1 ± 0.46] . (4.19) 

However, as for <jt(10°), small values of ag are required for n s < 0.6. If we use Eq.(4.19) together 
with Eq.(4.16) to constrain n s , the errors on the quadrupole are such that the range is not seriously 
restricted. (Again, we have ignored the asymmetry on the cosmic variance errors.) 

One can also use the correlation function data for given n s to determine the allowed range 
for ag. The correlation function (with quadrupole removed) and its cosmic variance are given by 
[50] 

C(0) = p e(cos 9) a\ t . (4.20a) 

and 

{[AC(9)] 2 ) = 2 £-Mp,(cos0) a 2 ,]' . (4.206) 

e>2 
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There are also correlations from angle to angle, so a matrix is more appropriate. As well, one 
should restrict the region of correlation function estimation to that actually used by the DMR 
team, which involved a cut in Galactic lattitude. This will increase the theoretical variance. In 
Figure 6, we compare our theoretical correlation functions, including their errors derived from 
Eq.(4.20), for the n s = 1 and n s = 0.4 cases with the DMR correlation function given in ref. 
[37]. We have fixed the amplitude of the theory curves by requiring that they give the DMR 
<7t(10°) = 1.09 x 10 -5 . If we vary this amplitude for fixed n s , then the theory will cease to agree 
with the data. Using the error bars that the DMR team give, and calculating \ 2 f° r the model 
fits to the data assuming the errors are independent and Gaussian (which they are not), we have 
constructed an allowed range for as which basically agrees with that derived from <tt(10°), but 
with slightly larger errors. A more precise treatment that takes into account the correlation in the 
variances of the theory C{9) and the influence of the extra correlation over pixel noise on the data 
C{9) error bars is needed to precisely pin down the allowed range. However, we are encouraged 
by the general agreement between limits derived from <7<r(10°), C(9) and the quadrupole. The 
DMR team derive the constraint n s = 1.1 ±0.5 from the correlation function data. Although it 
can be seen from Figure 6 that there is a slight preference for the n s = 1 case compared with the 
n s = 0.4 case, we do not consider that the n s = 0.4 case can be ruled out by this data alone. 

3) Large-scale Streaming Velocities 

There is another type of data that directly probes the amplitude of the mass density fluc- 
tuations as opposed to the fluctuations in galaxy or cluster number densities, namely large 
scale streaming velocities. From optical surveys, Bertschinger et al. [51] estimated the three- 
dimensional velocity dispersions of galaxies within spheres of radius 40h _1 Mpc and 60h _1 Mpc 
after the data had been smoothed with a Gaussian filter of 12h _1 Mpc, 

<r„(40) = 388 [1 ± 0.17] km s" 1 ; a v {m) = 327 [1 ± 0.25] km s" 1 , (4.21) 

which should be compared with the rms 3D velocity dispersions for power law CDM models (with 
errors calculated from the variance ([Acr 2 (40)] 2 )): 

a„(40) = 300 a 8 e im ^- n ^ [1±®>] km s" 1 ; a v (60) = 238 age 1 - 19 * 1 -"') [1±;||] km s" 1 . (4.22) 

The fits are good for < n s < 1. Although we do not regard these bulk flow estimates to be on 
as firm a foundation as the DMR measurement of <7t( 10°), it is interesting to note that the range 
suggested for as by the velocity data is similar, 

a s ~ 1.29e- 106(1 - ns) [l+f 5 ] , (4.23) 

provided n s is not very far from unity. It can be combined with Eq.(4.16) from <tt(10°) to yield 
a preferred value for n s of 1.07 (and a 8 = 1.4!), and a '2 sigma' lower bound of n s = 0.72. Using 
the 60h _1 Mpc c^-estimate gives a similar result. This constraint is so restrictive because the 
dramatic decrease in ag with decreasing n s from <tt(10°) more than offsets the increased velocity 
due to the enhanced large scale power. 

4) The Epoch of Structure Formation and Other Tests 

Given er 8 and the spectral index n s we can consider when structures of various types formed 
in the Universe. In Figure 7, we plot the range in linear rms density fluctuations a p (M) = 
((AM/M) 2 ) as a function of mass scale M allowed by Eq.(4.16). We actually calculate the rms 
fluctuations smoothed on a 'top hat' filtering scale Rth which is related to the mass by M « 
10 12 ' 4 (i?Ti?/h _1 Mpc) 3 . The range in Rth around R g = 0.5h _1 Mpc corresponds to the filtering 
appropriate for galaxy formation (top hat mass 10 n - 5 M Q ). The a p (M) shown are evaluated at 
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the current epoch if one extrapolates their growth by linear theory. This means that the rms 
fluctuations on the scale R g reach nonlinearity at a redshift somewhat above 

1 + z nl (R g ) = a p {R g ) « 6.2<r 8 e^ 1 "™*) « 7 . 2e -3-63(i-n s ) ^ ± a2] ? (424) 

where we have used Eq.(4.16) for <t 8 . Galaxies represent a much smaller fraction of space than 
that in typical fluctuations, but there is a lag between nonlinearity and complete collapse. These 
effects tend to cancel each other so Eq.(4.24) gives a first reasonable, although somewhat low, 
estimate of the redshift of galaxy formation. 

A better estimate of the redshift of galaxy formation is obtained in the following way. We 
take the observed luminosity function for galaxies [52] and assign an average mass-to-light ratio 
(M/L) for galaxies with luminosities above L. We then have, approximately, for the mass fraction 
in objects with luminosity greater than L, 



Q{> L) « 0.035exp(-L/L*)[(M/L)/(50h)]ft , 

where is a fitting parameter that gives the typical luminosity for a bright galaxy. The cor- 
responding mass is M = 6 x 10 11 h _1 [(M/L)/(50h)] L/L*. Therefore, the fraction of the mass 
in L* galaxies for the models we are considering is about a percent. Now consider the fraction 
of the mass in the Universe in collapsed objects with mass above 3 x 10 n M Q ; if we choose 50h 
for (M/L) and (M/L), this corresponds to the mass above L*/4, and the expression for f2(> L) 
above indicates that 2.7% of the mass should be in such objects. We thus determine the redshift 
at which the Press-Schechter mass function [53] for these models would predict that 2.7% of the 
mass in the Universe is in collapsed objects with mass above 3 x 10 11 M Q . The corresponding 
value for this redshift is just 30% higher than Eq.(4.24) and provides a better estimate of when 
pervasive galaxy formation would have occurred, 

(1 + z gf ) PS = 8.1(78 e-C""') « 9.5e- 3M3 ^- n ^ [1 ± 0.2] . (4.25) 

The power 3.63 is so large that even if we err on the conservative side by using Eq.(4.25) rather 
than Eq.(4.24) and take the upper limit, we obtain relatively strong limits on n s : 

n s > 0.63 , if z gf > 2 ; n s > 0.71 , if z gf > 3 . (4.26) 

A more careful analysis of star formation history would be required to improve upon these limits, 
but they illustrate that the amplitude factors allowed by the DMR data lead to strong limits on 
the spectral index to have galaxy formation occur early enough. Note that these bounds on n s 
are similar to those derived from the streaming velocities. 

A more powerful analysis of when objects of various masses form is provided by the hi- 
erarchical peaks method [54, 55], which identifies virialized potential wells with patches of the 
Universe centred on peaks of the density field that have undergone collapse, but solves the 'cloud- 
in-cloud problem' inherent in the original BBKS peak method [42] by merging small scale peak 
substructures into the dominant peaks that contain them. A mass function for dark matter halos 
at redshift z, n(M, z)dM, as well as detailed information about the spatial distribution of the 
halos, can be calculated. The objects found with this method have been shown to agree well 
with groups found in A-body calculations. Curiously, the mass function agrees reasonably well 
with that derived using the Press-Schechter approach [53], especially at the high mass end. This 
gives us some confidence in the validity of the Eq.(4.26), n s > 0.63, constraint. However, the 
Press-Schechter mass function has no strong theoretical justification [56] and cannot deal with 
the spatial distribution of objects. 

Since the total dark matter mass in galaxies is not directly measured, the mass function 
n(M) is of limited diagnostic use. On the other hand, the depth of galaxy and cluster potential 
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wells can be inferred from their internal velocity dispersion v. Therefore, in Figure 8 we show 
the number density of objects with velocity dispersion in excess of v, n(> v,z), for a variety of 
redshifts. The n s = 1 CDM model with as = 0.7 has roughly the right number of v = 200 km s _1 
halos at z > 3 to be a viable model of galaxy formation, and the number of clusters with 3D virial 
velocity above 1500 km s _1 roughly corresponds to the number of rich Abell clusters. Increasing 
<T8 for this model, as is suggested by the DMR data, might result in an excess of clusters with 
high velocity dispersions and thus high X-ray temperatures that may already be excluded by the 
X-ray data [55]. However, current indications from gravitational lensing observations in clusters 
[57] are that clusters exist with velocities in excess of v = 2000 km s _1 at z > 0.2, and a z ~ 0.2 
cluster observed with the X-ray satellite Ginga has an X-ray temperature of 13 keV [58], which 
translates into a»~ 2500km s _1 dispersion. It is also possible that cluster X-ray temperatures 
are below the values one would infer from the dark matter potential. Thus it may turn out that 
as ~ 1 will be preferred over 0.7 as the data improves. On the other hand, it is evident that 
cluster velocity dispersion estimates are easily contaminated by projection effects that always give 
overestimates [59], so the lack of v ~ 1500 km s _1 clusters in the as = 0.5, n s = 0.6 model cannot 
at present be used to exclude it. Thus, although it is universally agreed that the abundance of 
rich clusters as a function of velocity dispersion will be one of the most powerful measures of as, 
better data and extensive theoretical comparisons with the X-ray and optical data are required 
to test how strongly n s is constrained. The basic conclusion of the more complete analysis of 
ref.[55] is that, while one may argue that low amplitude models are not excluded by the velocity 
or temperature data, it seems quite unlikely that the errors in the X-ray flux and luminosity data, 
both for nearby and distant (z ~ 0.2) clusters, are so large as to allow these models to survive; 
explicitly, the n s = 0.6 CDM model with as < 0.5 is ruled out [55]. 

What even more strongly rules out the n s = 0.6 model, in agreement with the analytic 
argument constraining n s using z g j given above, is the lack of high redshift activity, in particular 
the paucity of halos with dispersion in excess of 200 km s _1 even as late as z = 2. These are the 
sites of bright galaxy formation. There are some interesting differences that appear at high z even 
with the modest change in slope from n s = 1 to 0.8, with as fixed: e.g., there would be an order of 
magnitude more v = 100 km s _1 'dwarf galaxies at z = 10 in the n s = 1 model than in the case of 
n s = 0.8. It has been argued [60] that only those dwarf galaxies with velocities above this number 
will survive the supernova explosions that occur when galaxies assemble themselves. Having some 
old cores of stable objects is probably a good thing rather than a bad thing, since they could be 
the birthplaces of quasars, but because of uncertainties in modelling the gas dynamical behaviour 
of forming galaxies and of the intergalactic medium one cannot be sufficiently definitive about 
the high z consequences of a theory to select one model over the other at this stage. 

Another test which has been used to argue that as ^ 0.6 and which therefore favours n s < 1 
models is the velocity dispersion of pairs of galaxies over separations of order a Mpc [61]. In the 
early A r -body simulations of n s = 1, as = 1 CDM models, the pair velocity dispersion of dark 
matter halos on these scales was found to be much higher than the velocities of galaxies inferred 
from redshift surveys. However, Carlberg and Couchman [62] computed an n s = 1 CDM model 
in which the relative velocity of galaxies was much less than that for the dark matter, an effect 
termed 'velocity bias'. Coincidently, they chose as = 1.17, the value suggested by DMR. Although 
how effective this velocity biasing can have been at lowering the pair velocities is a matter of much 
debate, smaller n s will obviously help to ease the problem. 

Experimental upper limits on small and intermediate angle anisotropies in the microwave 
background can also be used to constrain the index n s , but require detailed computations along 
the lines of those given in ref. [50] and we shall not undertake them here. We note however that 
the pre-COBE limits on anisotropy were already strong enough to place constraints of n s > 0.6 
for <jg = 1 and n s > 0.3 for a s = 0.5 [47] at the 90% confidence level, and the constraints from 
an earlier DMR limit [63] also gave similar values. (For other previous discussions of power law 
CDM spectra, see [38,39,67].) 
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5) The Role of Gravitational Wave Modes 

Stimulated by the DMR results, other groups have been independently considering inflation- 
inspired power law spectra [68, 69]. Davis el al. [69] have pointed out that, although gravitational 
wave modes are generally small for nearly scale invariant spectra [70], for n s Cl this conclusion 
may not hold, amplifying upon the work of Abbott and Wise [71]. Although gravitational waves 
do not make an important contribution for natural inflation, they are significant for power law 
and extended inflation models. 

We first sketch why they can be ignored in natural inflation. During inflation, the same 
zero point quantum fluctuation phenomenon which leads to the inflaton density perturbations 
also leads to statistically independent gravitational wave perturbations. If h + and h x are the 
two linear gravitational wave perturbations, then </?+. x = Mpi h +jX /\/16tt behave just like single 
massless scalar field degrees of freedom as far as fluctuation generation is concerned. Each of 

1/2 

the fields v?+ iX of comoving wavenumber k have power spectra P<p + x (k) equal to the Hawking 
temperature H/(2ir) when k = Ha, just as the inflaton fluctuations do, except that they are not 
amplified during subsequent evolution. With the factor given above, we therefore have for the 
total gravitational wave power, 

P GW = i p h+ + Ph x ] 1/2 = V^M pi - 1 H/{2t:). The ratio of the 
gravitational wave power spectrum to adiabatic metric perturbations, as encoded in the spectrum 
Pf , at horizon crossing is therefore 

p g 'w = vWie^l 

pV2 3M Pl H ' 

where the y/2 comes from the 2 independent GW polarizations that can be generated. Using the 
WKB values at horizon crossing usually gives accurate estimates of final fluctuation amplitude 
[11]. For natural inflation, and using the slow roll approximation and the results of Sec. IV. 1, we 
have 



P g'w _ 2^2 _ 2^2 / M Pl 2 x 1/2 W M Pl 2 \ 



2 



M P i 



1 



-1/2 



(4.27) 

Thus the gravity waves are exponentially suppressed relative to the adiabatic scalar fluctuations 
of the inflaton over the observable large scale structure waveband. In particular, for / < M P i, 
this ratio is less than 0.04 for modes with wavelength equal to the current Hubble radius. On the 
other hand, for power law inflation with an exponential potential, the ratio is 



Pj/w _ 2^2 _ 2^2 
P V2 3 



1 + 



n. 



-1/2 



(4.28) 



which can be quite favourable to the tensor modes if n s is sufficiently small. 

The amplitude of gravitational wave modes decays by directional dispersions as the modes 
re-enter the horizon, just as waves in any relativistic collisionless matter do [43], whereas the 
adiabatic fluctuations maintain a constant gravitational potential. Before the gravitational wave 
structure disperses however, it influences the microwave background through the Sachs- Wolfe 
effect. A number of authors have calculated the magnitude of this effect [70,71]. We denote 
the ratio of tensor to scalar contributions to the radiation field multipole moments o,lm by Al. 
Abbott and Wise [71] show that this ratio is not very sensitive to the multipole moment L. Davis 
et al. [69] use the results of [70,71] to get the ratio for the quadrupole value; in our language, this 
is 

^ - - 3.9§^ . (4.29) 

°T i=2 P c 
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To estimate the correction for power law inflation, we shall assume Al is A2, which, using 
Eq.(4.28), is therefore 



A 2 « 3.7 



1 + 



1 - n s 



-1/2 



(4.30) 



The value of ctt( 10°) given in Eq.(4.14) should be multiplied by [1 + A 2 ,] 1 / 2 . Thus, the range 
for o"8 as a function of n s is lowered substantially as a result of the inclusion of gravity waves, 
as we have shown in Figure 5; e.g., a s drops by a factor of 1.8 for n s = 0.6. This makes the 
already strong constraints we have derived significantly stronger. The n s -constraint we derived 
by requiring that galaxies form early enough in the theory, n s > 0.63 for z g f > 2, changes to 
n s > 0.76 for power law inflation; similarly, the bound n s > 0.71 from the requirement z g f > 3 
now becomes n s > 0.82. Also, the '2 sigma' streaming velocity limit of n s > 0.72 increases to 
n s > 0.89. 

For the chaotic inflation potentials used above, we have 



pl/2 
GW 
pl/2 



N T (k) 



-1/2 



A 2 « 2.63^ N T (k) + 



■1/2 



(4.31) 



hence gravity waves diminish 8 by only 11% for a (p 4 potential, and by 5.5% for a 4> 2 potential. 
Slightly higher values are obtained if we use a power law inflation formula with n s = 0.95 and 0.97, 
respectively. Again motivated by COBE, various authors have been looking at the gravitational 
wave contribution in these conventional inflation models anew [69,72]. 

It is clear from this discussion that if one could unearth the gravity wave component of 
anisotropy from the adiabatic component, it would 8not only allow a strong discrimination among 
models, but it would also rule out natural inflation, which predicts no component whatsoever. 



6) Discussion 

Since our n s > 0.6 limit comes from a variety of arguments, we believe it is quite robust. Thus, 
unless the errors in the analysis of the large scale clustering observations are larger than currently 
estimated, a fluctuation spectrum with broken scale invariance that has a slowly changing spectral 
index over the range from ~ 10 — 10 4 Mpc cannot be the solution to the extra power dilemma 
that the CDM model faces. However, the allowed values of n s > 0.7 can help to ease the 
requirements on some of the extra power fixes proposed in the literature (e.g., [45]). 

Motivated by the DMR results and the many prospects for broken scale invariance in inflation, 
Cen et al. [68] have very recently undertaken combined hydrodynamical and iV-body calculations 
of CDM models with n s = 0.7 and have independently come to a number of the conclusions we 
have about such models, namely that they help but do not fully solve the large scale structure 
dilemma. 

Finally, our limit on n s can be translated into constraints on the parameters of inflation 
models that give rise to power law spectra. For example, it gives a very strong constraint on 
the effective value of uj, the Brans-Dicke parameter which arises in extended inflation models. 
When the effect of tensor waves is included, the z g f > 2 constraint, n s > 0.76, becomes uj > 17, 
near the upper limit of the range uj < 25 required for successful inflation [73] in most versions 
of this theory. Indeed, a closer examination [74] of the upper bound on uj, which arises from the 
requirement that large bubbles do not produce an excessive microwave anisotropy, suggests that 
in fact uj < 18 is required if the dark matter is cold. (This number might even be slightly lower, 
since it is based on the older COBE data.) Combined with our lower bound on ui, this limit would 
leave little room for most extended inflation scenarios. For natural inflation, from Eqn.(4.2), the 
constraint n s > 0.63 translates into a lower bound for / of 0.33Mp;. This is comparable to the 
constraint (3.37) from reheating. 
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V. CONCLUSIONS 



We have studied an inflation scenario inspired by particle physics models with weakly self- 
coupled (pseudo-)scalars such as the axion. With the requisite mass scales, which can emerge 
dynamically for plausible choices of gauge groups, PNGB inflation appears to be robust in the 
sense that it arises in the simplest class of models, with a potential of the form (1.1). We have 
shown how these models can arise in a variety of theoretical settings, and indeed that superstring 
models already in the literature come very close to providing the desired mass parameters for 
natural inflation. Although the tendency of higher dimension operators on PNGBs arising from 
wormhole effects, for example, would be to increase A, we discussed quite plausible ways in which 
the upward movement can be exponentially suppressed, so our model retains its naturalness. 

We numerically and analytically studied the cosmological dynamics of the inflaton field, and 
derived several constraints on the two-dimensional parameter space (/, A). The allowed band of 
parameter space includes models which have more relative fluctuation power on large lengthscales 
than the standard scale-invariant spectrum. We have studied in depth the consequences of these 
power law initial fluctuation spectra for large-scale structure and the microwave background 
anisotropy. We find that models with n s < 0.6 are required to fit the large-scale galaxy angular 
correlation function w gg {9) observed in the APM survey, but the recent COBE results require 
a rather small amplitude for these models to be consistent, as 0.5. This makes the epoch of 
galaxy formation uncomfortably recent and predicts large-scale flows of relatively small amplitude. 
Turning this argument around, combining the COBE results with the requirement of sufficiently 
early galaxy formation and large-scale flows of the inferred amplitude leads to the constraint 
n s > 0.6. For natural inflation, this implies / > 0.3Mpi, virtually the same bound as we get from 
the reheating constraint. Although the simple expedient of reducing n s does not, by itself, solve 
all the large scale structure dilemmas for the CDM model, it can be combined with other ways 
to explain the extra large scale power [45], for example, by introducing into the CDM model a 
neutrino with a mass of a few eV, a nonzero cosmological constant (Mpi A/87rh = 0.2 with CDM 
fits for n s = 1), a smaller Hubble constant (h ~ 0.4), a larger baryon abundance, or by simply 
supposing that galaxies are distributed on large scales somewhat differently than the mass so 
that the linear biasing assumption of Eq.(4.8) is invalid. We conclude that inflation with pseudo- 
Nambu-Goldstone bosons offers an attractive model for generating curvature fluctuations whose 
gravitational instability can lead to all of the cosmological structure we observe around us, even 
if the spectrum is nearly scale invariant. 
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APPENDIX: APPROXIMATION OF INTEGRALS 



In this Appendix, we demonstrate the validity of the results presented in §111. C on the a 
posteriori probability of inflation. The main difficulty is that the simple logarithmic form (Eq. 
3.11) for the number of e-foldings as a function of yi = <fii/ f does not hold for large y\ [i.e., for 
2/i > !)• We should thus write the integral / (see Eq. 3.31) in the form 



. (Al) 

' ' 1P 3JV( V1 ) 

\\B -.vi- 
lli §111. C, we argued that when 3B > 1, the integral can be approximated by the first term above, 



3B 

We now calculate the relative error suffered in making this approximation. We first note that the 
number of e-foldings N(yi) is a strictly decreasing function of the starting value y\. In particular, 

N( yi )<N(l) = A Vj/i€[l,4 (A3) 

8We thus obtain a bound on the second integral in Eq. [Al], 



£ d yi e 3N ^ <e 3A (ir-l). 



(A4) 



This contribution to the error is always positive, whereas the other contribution [namely — e 3A / (3B- 
1)] is always negative. The total error E is therefore bounded from above by 



E< e 3A 



7T- 1 - 



35-1 



(A5) 



The total error is also bounded from below by the second (negative) term alone, so we obtain the 
relation 

-1 < E(3B - l)e~ 3A < 3B(ir - 1) - tt, (A6) 
and hence the relative error £ = E/I is bounded by 

£ < e 3B ~ l x max{l, 3S(vr - 1) - vr}. (A7) 

This error is always sufficiently small for the cases of interest. For example, for / w Mpi, 
3B ~ 48vr(//Mp2) 2 « 48vr, then e < y? ax ~ 0.6, and hence £ < 2 x 1CT 31 . For the other end 
of the mass range of interest (i.e., for / near f c = 0.06Mp;), let 3B — 1 = 5 where 5 is a small 
positive number. In this regime y™ ax ~ 10 -60 and hence £ < 10 -605 . The error is thus completely 
negligible until 5 becomes smaller than 1/60 or so, that is, until / is very close to f c . 
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FIGURE CAPTIONS 



Figure 1: Plot of various field values and parameters vs. / /Mpi. The upper curves show that 
our estimate fa/ f (eq.3.5) for the value of the field at the end of the SR epoch (dotted) is very 
close to our numerical result <p e nd/ f for when inflation ends (solid); the middle (dotted) curve 
shows log(4>™ ax / f) , the largest initial value of the field consistent with 60 e-folds of inflation; the 
lower (dashed) curve shows the density perturbation constraint (3.23) on the scale A [plotted as 
log(A/Mp/)], assuming the bias parameter b g = 1. 

Figure 2: Results of the numerical integration of the scalar and gravitational equations of motion. 
The number of inflation e-folds N(fa) of the scale factor is shown as a function of the initial value 
of the scalar field, fa, for different values of the fundamental mass scale, f/Mpi = 0.05, 0.07, 0.1, 
0.2, and 0.5. 

Figure 3(a): Power spectra for CDM models with variable spectral indices n s are plotted against 
comoving wavenumber k (referred to current length units). There is progressively more large scale 
power as n s decreases through the values n s = 1, 0.6, 0.4, 0, -0.4, and -1 shown in the figure. The 
lines under the labels (whose vertical placements are arbitrary) indicate approximate regions in k 
space that various probes of structure are sensitive to, such as: microwave background anisotropy 
experiments of large angle (COBE, [37]) and of intermediate angle (e.g., SPole is a 1° experiment 
[64]); clustering observations for galaxies in the APM Galaxy Survey (w gg , [46]) and the QDOT 
redshift survey [65] and for clusters (£ cc ) [48]; and large scale streaming velocities (LSSV) [51]. 
The hatched region denotes the region that the power spectrum must pass through to explain the 
APM angular galaxy correlation function data [46]. The spectra are in units of a$. 

Figure 3(b): Galaxy power spectra derived assuming linear dynamics (appropriate for k~ l > 
5<78h -1 Mpc and large scale linear biasing), in units of b g ag, for CDM models with variable 
spectral indices n s . This region of the spectrum is highlighted because that is where the large 
scale structure data exists. The hatched region is the APM region of (a), while the points denote 
the power spectra estimated from the QDOT redshift survey [65] and the IRAS 1.2 Jansky survey 
[66]. The biasing factors for the (slightly) different types of galaxies probed by the APM, QDOT 
and 1.2 Jy surveys could explain the differences in these results. There are indications that 
b g = OSa^ 1 is needed for the 1.2 Jy galaxies [66], while b g = a^ 1 describes the APM survey well, 
and this relative factor is enough to bring the required power spectra into line; i.e., the n s = 
0.2-0.6 range is also preferred by the 1.2Jy data if b g = OSa^ 1 , while the n s = 1 curve falls below 
the data error bars. 

Figure 4: The models of Fig. 3(b) (with n s = 1, 0.8, 0.6, 0.4,..., -1) are compared with the 
angular correlation functions determined from the APM Galaxy Survey [46] scaled to the depth 
of the Lick catalogue, at which 1° corresponds to a physical scale of ~ 5/i _1 Mpc (dots). No 
nonlinear corrections were applied to the theoretical power spectra, but for angular scales above 
~ 1° and for amplitude factors a 8 < 1, the linear approximation is accurate [44]. The theoretical 
curves are in units of (b g ag) 2 . The straight line gives the angular correlation that would result if 
the behavior of the spatial correlation function observed over distances r < 10/i -1 Mpc, £ ~ r -18 , 
were extended to large separations. The hatched region corresponds to the allowed region once 
corrections for systematic errors are included. The data therefore suggests < n s < 0.6 is needed 
for the CDM model if biasing is linear on large scales. 

Figure 5 The range (with '1 sigma error bars') of the amplitude parameter ug for a standard 
CDM model in the limit that Vtp = as a function of the power law slope n s , using the constraint 
from the rms fluctuations on 10° in COBE's DMR experiment. Both the theoretical variance and 
the quoted experimental error are included in the error bars, which are in total about ±20%. The 
values of as drop by a further ~ 15% when tip ~ 0.06 is used rather than the zero used here. If 
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the correlation function data of Fig. 6 is used to determine the amplitudes, a similar constraint 
curve arises. Also indicated is the range in n s suggested by the APM angular correlation function 
data, the '1 sigma error bars' on n s derived using the correlation function by the DMR team, and 
the as range (dotted lines) that encompasses the values that have been advocated for the n s = 1 
CDM model by different workers, e.g., as = 0.4, 0.55, 0.65, and 1.2 in [61,44,54,62] respectively. 
The dashed curves give the allowed range for o"g(n s ) when gravitational wave modes are included 
for power law (and extended) inflation. For natural inflation, the deviation from the solid curves 
is infinitesimal. 

Figure 6: Comparison of the DMR 53A + B x 90A + B cross correlation function with the 
quadrupole removed [37] with the theoretical predictions, including variance, for (a) n s = 1 and 
(b) n s = 0.4 spectra. The as amplitudes shown have been set by requiring the angular power 
spectrum to reproduce the rms fluctuations on 10°. Clearly, although the data is somewhat better 
fit by the n s = 1 rather than the n s = 0.4 model, one cannot strongly distinguish between the 
2 models on the basis of shape alone. The experimental errors plus theoretical variance in the 
quadrupole amplitude are sufficiently large that one cannot use the comparison of the quadrupole 
with <tt(10°) to effectively constrain n s . The strongest restriction comes from the consequences 
of low as for n s < 0.6 for structure formation. 

Figure 7: The linear rms fluctuations averaged over spherical regions of radius Rth are plotted 
as a function of the mass M « 10 12A (R T h/ h _1 Mpc) 3 M , for CDM models with n s = 0.4, 0.6, 
0.8, and 1 (with n s increasing as one moves vertically up the figure). The error bars show the 1 
sigma range in spectrum normalization as a result of DMR and cosmic variance errors in or (10°). 
Although Fig. 3(a) shows that n s < l spectra have more power on large scales and less on small 
scales than n s = 1 models with the same as, when as{n s ) determined from COBE is used, the 
amplitude for n s < 1 is less on all mass scales. The extreme problems with the n s = 0.4 model 
and the marginality of the n s = 0.6 model are evident from this graph alone. 

Figure 8: In (a), we show the number density of collapsed objects with 3D virial velocity in 
excess of v for the CDM model with spectral index n s = 0.8 and for the value of the amplitude 
parameter as = 0.7 (indicated by the DMR 0-^(10°) data for this model). The densities are 
shown as a function of redshift z, with z decreasing as one moves to the right in the figure. The 
velocities in the hierarchical peaks method [55] used for this computation could be larger by an 
amount given by the error bar labelled by 'v range'; these error bars are explicitly put on the 
z = curve. The number densities shown should be compared with the abundances indicated 
by the horizontal lines and velocity dispersions indicated by the downward arrows: for 'bright' 
galaxies, ~ 10~ 2 ( h _1 Mpc)~ 3 with v ~ 220km s -1 , for rich clusters, ~ 6 x 10~ 6 ( h _1 Mpc)~ 3 with 
v ~ 1500km s _1 , and for at least one object between us and redshift 2, ~ 10 _9 ( h _1 Mpc) -3 with 
v ~ 2500km s _1 , according to the Ginga X-ray satellite team [58]. In (b), we choose the DMR 
1 sigma upper bound as = 0.5 for n s = 0.6; even so, the number of 'bright galaxy' halos is too 
small by z = 2. In (c), we plot the densities for n s = 1, using the DMR 2 sigma lower bound 
a 8 = 0.7 for the amplitude, to facilitate comparison with (a). The number densities of model (c) 
accord reasonably well with the hierarchy of objects in the Universe. There is little to distinguish 
between the n s = 1 and n s = 0.8 models with the same erg. To explicitly show this, we also 
plot with light solid curves the tails of the z = abundances for cases (a) and (b). The third 
light curve, also for z = (the highest curve at large v), shows the effect of increasing as to 1 
for the n s = 1 model, closer to the number indicated by DMR. Although this may lead to too 
many clusters with higher X-ray temperatures than observed [55], ug = 1 does help to explain 
the Ginga event. 
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